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We introduce Wilson's, or Polchinski's, exact renormalization group, and review 
the Local Potential Approximation as applied to scalar field theory. Focusing on 
the Polchinski flow equation, standard methods arc investigated, and by choos- 
ing restrictions to some sub-manifold of coupling constant space we arrive at a 
very promising variational approximation method. Within the Local Potential 
Approximation, we construct a function, C, of the coupling constants; it has the 
property that (for unitary theories) it decreases monotonically along flows and 
is stationary only at flxcd points - where it 'counts degrees of freedom', i.e. is 
extensive, counting one for each Gaussian scalar. 

In the latter part of the thesis, the Local Potential Approximation is used to 
derive a non-trivial Polchinski flow equation to include Fermi flelds. Our flow 
equation does not support chirally invariant solutions and does not reproduce 
the features associated with the corresponding invariant theories. We solve both 
for a flnite number of components, N, and within the large A'^ limit. The Legendre 
flow equation provides a comparison with exact results in the large A^ limit. In 
this limit, it is solved to yield both chirally invariant and non-invariant solutions. 
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Chapter 1 



Introduction 



1.1 Motivation and Structure 



The purist aspires to a reliable method of performing non-perturbative calcula- 
tions in quantum field theory whereas the practically minded physicist desires an 
efficient cost effective procedure of accurately generating physical quantities. We 
will argue that the methods presented here offer us hope of attaining these desires 
and aspirations. The strength of the scheme discussed exhibits itself in two impor- 
tant features. First the method can be improved in a systematic way producing 
a sequence of converging approximations and second, there is no small parameter 
to control the approximation (unlike perturbation theory), allowing us to reach 
and test new and interesting areas of physics. In particular we demonstrate that 
the method produces competitive results in the context of scalar field theory, 
both numerically and theoretically via critical exponents and Zamalodchikov's 
C-theorem respectively. Many authors use these methods for alternative reasons, 
e.g. to calculate effective potentials. The successes presented here motivate our 
attempts to extend the formalism to include Fermions. 
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We begin the thesis by introducing the essential concepts of the renormahzation 
group methods from a field theory perspective. The inevitable parallels with 
theories of critical phenomena are drawn with the usual lengthy and detailed dis- 
cussions left to reference [1-5] . Then the methods applied to scalar field theory are 
reviewed and compared, paying particular attention to a Z2-invariant one com- 
ponent field theory in three dimensions. We develop a new variational approach 
which is applied to the scalar field theory under consideration and compared to 
the more standard techniques. We give a brief introduction to conformal field 
theory and review Zamalodchikov's C-theorem. Then, completing the defence of 
the Local Potential Approximation, we show how an appropriate generalisation 
of Zamalodchikov's C-function can be found within this scheme. The thesis fin- 
ishes with an outline of progress made and problems associated with attempts to 
extend the work beyond purely scalar field theory. In particular we compare the 
Local Potential Approximation in the context of Fermionic field theories with ex- 
act results obtained in the large N limit. This leads us to an understanding of the 
corresponding results at finite N. We will work in Euclidean space throughout, 
unless stated otherwise. 



1.2 Renormalization 



Prior to the main discussion of the thesis we wish to stress the modern point of 
view of renormalization as distinct to that presented in older textbooks. Here, 
we define renormalization as the procedure of re-expressing the parameters which 
define a problem in terms of some other, perhaps simpler set, while keeping 
unchanged those physical aspects of the problem which are of interest [4] . 
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1.3 The Ising model 



A classical theory of spin-spin coupling in a ferromagnet can be described by the 
Hamiltonian, 

H{s,h) = -^Vn,mSn-Sm + h^Sn (l-l) 

n,m n 

where h is the external magnetic field and Vn,m is the coupling of spins s„ and 
Sm at sites x„ and Xm respectively. Neighbouring sites in the lattice of spins are 
separated by a distance a. If we assume the lattice is /^-dimensional hyper-cubic 
and we only have nearest neighbour interactions, we can write 

-J2^n,mSn.Sm = Kj2[J2'^Sn+f,~S^)^-2Dsl] (1.2) 

n,m n \ fj, / 

where /i = 1, - ■ ■ ,D. Within the Ising model s is a scalar which takes the values 
s = ±1. An effective]^ description of the partition function is given by 

Z{h,f3) = I n^w(sn)e-^^(^''^) (1.3) 

where (3 = -j^^ and p(s„) is a weight factor describing the averaged local, micro- 
scopic properties of the spin. A choice like 

p{Sn) OC 6-'=^"-^^" (1.4) 

gives a convenient description. The partition function (1.3) may now be written 

as 

Z{K, u,X,h) = JY[ ds^e-^^'^'^'"'^''^) (1.5) 



^we replace the discrete spins by continuous variables which by universality will produce the 
critical behaviour of the Ising model 
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Figure 1.1: Symmetry breaking 



where the effective Hamiltonian is given by 



H{s; K, u,X,h) = J2( Y.{s^+, - s^f + i^{P)sl + A(/3)4 + hs^ (1.6) 



and v{(5) = c(/3) - 2D. Consider the potential v{(5)s^ + \{(5)s^ (figure 1.1). A 
ferromagnetic transition (symmetry breaking) occurs at y{(3c) — 0. The value of 
T, for which v{l3c) = 0, is called the critical temperature, Tc- 

Minimising the effective Hamiltonian in (1.6) with i/{f3) > leads to a ground 

state where all the s„ = 0, leading to zero magnetisation, 



However the ground state for i/(/3) < corresponds to all s„ aligned with s„ = 



n 




(1.7) 



n 
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, leading to magnetisation, 




The Ginzburg Landau theory of ferromagnetic transitions assumes that K{P), 
and A(/3) are smooth functions of temperature. Thus we might write 

cx (/?-/?,) (1.9) 

for (3 ^ j3ci leading to 

< s >oc ^/(3-/3c (1.10) 

for (3 > /?c. This shows the typical non-analytic behaviour at a phase transition, as 
indicated &i (3 — (3c in figure 1.2. Near the critical point, translational invariance 
implies that the correlation function will behave as 

< (S„- < S >){Sm- < S >) >~ \Xn - T^-^"^"^"^ (l-ll) 

for the range a <ti \xn — Xm\ ^ ^{h,l3). The lack of any dependence on a 
fundamental length scale reflects that at criticality the system is scale invariant 
(to be discussed in greater depth later). The correlation length, ^, can be thought 
of as the size to which we can reduce the system without qualitatively changing 
its properties. The anomalous dimension, 77, is effectively defined by (1.11). We 
will discuss this quantity in greater detail later. 
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Figure 1.2: Magnetisation as a function of temperature for the Ising model with 
zero external magnetic field 

1.4 Lattice Field theory 

Consider the Lagrangian for one component 0*^ scalar field theory in D-dimensional 
Euclidean space, 

+ ^mV + ^0^ (1.12) 

We choose to regularise the field theory by discretising spacetime using a finite 
volume hypercubic lattice. This leads to a direct comparison with the Ising model 
of a ferromagnet. We let 

X — > a;„ = ane^, (1-13) 

(f){x) ^ 0n = Hxn), (1-14) 

1 1 

Qj a 

and the measure 

V 

Vct>^Y[d(t>n (1.16) 

n 

where denotes D orthonormal vectors and a is the lattice spacing. 
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By scaling the fields, sources and coupling constants, such that they are dimen- 
sionless, 

A = /a^-4, 0' = ^a5(^-2)0, J' = (^a^(^+2)J, (1.17) 
the discretised action can be written 

5(0', J') =Y.{\ - ^'nf + \^'<^'<^n + ^01' + ■ (1-18) 

If we compare this with the Hamiltonian for the spin system (1.6) we see that at 
the classical (or Ginzburg Landau) level the continuum limit of the dimensionless 
lattice field theory (a 0) corresponds to approaching the critical point of the 
ferromagnetic transition (z/(/5) 0) because of the identification 

mV ~ v{f5). (1.19) 

Furthermore, we associate the even symmetry of the lattice field theory (0 — >■ 
—0)0 with a 'fiipped spin' symmetry of the spin system (s — > — s). Thus the 
continuum field theory can be viewed as a critical classical spin system. 



1.5 Universality 

The association between the lattice 0"^ theory and the Ising model may appear 
contrived, however it is found that different physical systems form 'universality 
classes', in the sense that quantum field theories can be used to describe critical 
phenomena for a range of models. Although there is no proof, there is consid- 
erable evidence that 0{N) symmetric scalar field theories correspond to various 
universality classes according to the value of N [1] ; 

^later we will be interested in a Z2-invariant one component scalar field theory 
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N—0: Critical behaviour of polymers. This is only formally defined in the limit 
— > as noted by Gennes [6] . 

N—1: Liquid-vapour transition and alloy order- disorder transition. Observe that 
we have no internal symmetries and the liquid- vapour transition can be modelled 
by particles living on a lattice, allowing the occupation of each site to be either 
zero or unity (equivalent to the Ising model) [7]. 

N=2: He^ superfluid phase transition and planar ferromagnets. The planar ferro- 
magnet has symmetry in a plane corresponding to an 0(2) invariant scalar field 
theory. 

N—3: Ferromagnetic phase transitions. Here we observe a three dimensional 
symmetry corresponding to an 0(3) invariant field theory. 

N—4: Chiral phase transition. It has been postulated to correspond to the chiral 
phase transition for two fiavours of quarks. 

Arguments presented in the following analysis provide a phenomenological expla- 
nation of universality. 

1.6 Kadanoff blocking 

We now drop the distinction between the Ising model and the lattice field theory 
and use a renormalization procedure known as Kadanoff blocking [8] to reduce 
the number of degrees of freedom. Our toy model Hamiltonian will be 

^ = E f E \^^n^. - ^n? + + Ml] . (1-20) 
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We begin by dividing the original lattice into blocks of size where s is an 
integer, and define an average field in each block Bs{n') labelled n' (figure 1.3) 

by 

i E <t>n. (1.21) 



ntBs (n') 



We end by scaling the blocks back to the original size 

Xs = -^ (1.22) 

and 

0, = (1.23) 

to ensure that the short distance properties of the correlation functions (physics) 

are left unchanged. This can be seen by considering the effect of the block 
transformation on the correlation functions: 



n'W/ - ^< I] I] '^m> (1.24) 

neBs(n') rraeBs(O') 
! < 0,„/0o > (1.25) 



since all the correlation functions in (1.24) are at essentially the same distance 
if n' is very large. Then, from (1.11) we see that 

< <t^n<t>o (1.26) 



and thus for the correlation functions to be left unchanged by blocking we must 
simultaneously scale the fields as indicated by (1.23). 

The new Hamiltonian, ifs(0s), will not be identical to the original, (1.20), but 
will contain infinitely many additional arbitrarily complicated terms including 
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Figure 1.3: Kadanoff Blocking 

for example 

(^n+M ~ ^^n)^0n> ^n+M^n ^nd - 20„ + (1.27) 

Thus, given that we will consider successive blockings, it is natural to start with 
a completely general action including all terms which obey the symmetries of the 
theory under consideration. In general these terms will be infinite in number. 
Then, each blocking can be viewed as a mapping of the 'coupling constant space' 
onto itself, called the renormalization group transformation, denoted T{s). We 
can see the group structure exphcity through T{si)T{s2) — T(siS2), although 
strictly speaking we only have a semi-group, since there will not necessarily be 
an inverse transformation. 

In terms of the ferromagnet we imagine looking at a sample through a micro- 
scope such that our eyes can see spin variations down to a certain size. Then 
Kadanoff blocking represents the operation of decreasing the magnification of the 
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microscope by the factor s. 

1.7 Effective Lagrangians 

The lattice we discussed earher is simply a method of regularisation which cor- 
responds (more or less) to introducing an overall momentum cutoff (Aq ~ ^). 
Thus, Kadanoff blocking can be interpreted as integrating out high momentum 
modes0 and scaling the cutoff back to its original size, A = Aq, such that the 
physics remains unchanged. However, it is simpler and equivalent to ensure that 
all variables are 'measured' in units of A, i.e. we change variables to ones that are 
dimensionless, by dividing by A raised to the power of their scaling dimensions. 
Then the action of lowering A reproduces the scaling steps in (1.22) and (1.23). 

Here we are considering directly a continuum effective Lagrangian describing 
physics over a limited range of energies, i.e. momentum less than Ag. The most 
famous example of an effective theory is the Standard Model which produces good 
results below the scale of Grand Unification, where new physics occurs. Through- 
out this thesis we will be dealing primarily with the Legendre and Wilsonian (or 
Polchinski) effective actions. Begin by considering a general effective action, S'^-^^ , 
which we assume to be valid below the overall momentum cutoff, Aq. Below Aq 
the effective Lagrangian will be described by an infinite set of couplings. Then, to 
look at physics at some energy scale E' ^ A^ we lower the cutoff to a scale A ~ 
and allow the couplings to flow in order to keep the physics constant (partition 
function invariant with respect to A). Then the action will evolve according to a 
flow equation of the general form 



A 




dA 



■'introducing a secondary momentum cutoff, A 
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Conventionally we parametrise flow equations in terms of the 'renormalization 
time', defined by 

t = ln(^] (1.29) 



so that t runs from to oo as A runs from Aq to 0. Generally we will be interested 
in solving in the continuum limit, corresponding to Aq — > cxd, where we require 
the solutions of interest to be insensitive to Aq. 



The effective theory is said to be renormalizable if we can calculate physical 
processes once we have determined a finite number of couplings, known as rele- 
vant. Relevant couplings , denoted g^, flow away from some 'fixed point' value, 
g^. The other, irrelevant couplings, g^, flow toward some 'fixed point' value, g^. 
The renormalization (semi-) group moves the action along the trajectories in the 
coupling space. 



1.8 Fixed points and Perturbations 

If the bare action (at A = Aq) is chosen in the neighbourhood of the fixed point 
in such a way that the relevant couplings are their fixed point values {g^ = g^) 
and the irrelevant couplings are close to their fixed point values then the action 
will flow into that fixed point, -S';^-''-^, defined by 

- f[Sl"] = 0. (1.30) 

A fixed point is associated with a critical surface in the coupling space which 
contains all the bare actions which flow into that fixed point. Figure 1.4 illustrates 
the neighbourhood of a fixed point with the critical surface and renormalization 
group trajectory associated with it. In practice this space is infinite dimensional. 
Similarly, if we begin with the ferromagnet at its critical point, T — Tc, and 
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decrease the magnification sufficiently by Kadanoff blocking, eventually we shall 
observe no further change {i.e. we reach a fixed point). 

We enquire about the stability of the fixed points by perturbing about the fixed 
point solution. For example, later we will expand a dimensionless potential (for 
dimensionless scalar fields) as 

V{cf>,t) = K(0) + e^/„(0)e^"* (1.31) 

n 

where K is the fixed point solution under consideration and the perturbation is 
considered for small e. We immediately notice that A„ > are associated with 
relevant directions in the coupling space, in which the perturbations grow with t 
and similarly A„ < are associated with irrelevant directions. Fixed points are 
characterised by the A„ associated with them. If = we can to first order 
add any amount of the corresponding 'marginal perturbation' to the fixed point 
without introducing any t dependence. Such a scenario can result in lines of fixed 
points, and does in a few specific cases to be discussed. 

Substituting (1.29) into (1.31) we find that the dimensionless perturbation around 
a fixed point is given by 

SV^eY,fn{<P)(^T . (1.32) 

n ^ ^ 

However, more generally the fiow equation (1.28) implies that we can write 

A J^(^Fa = -L{SVa) (1.33) 

where L is a linear operator acting on 6Va. Thus from (1.32) we see that the A„ 
form a discrete set of eigenvalues for L, where the /„(0) are the corresponding 
eigenfunctions. 
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Figure 1.4: A two dimensional projection of the flow near a fixed point 

Since we exchange all dimensionful variables for dimensionless ones using A, in- 
dependence of A would imply lack of dependence on any scale at all. Thus a 
fixed point, defined by (1.30), corresponds to a scale invariant theory and thus is 
associated with either divergent or vanishing correlation lengths. Hence a fixed 
point is massless (or infinitely massive) respectively. A massive theory leads to 
no propagation and is thus considered less interesting. Later, we will discover a 
trivial fixed point known as the Gaussian which represents a massless quantum 
field theory. The value of the coefficient of 0^ will generally be non-zero at a mass- 
less fixed point (with the exception of the Gaussian), but will be precisely the 
required 'counter-term' to cancel the remaining quantum corrections, leaving the 
theory overall massless. Furthermore, the massless fixed point must correspond 
to the continuum hmit (Aq oo), for otherwise A^ would set the scale. 

Now we are in a position to make the concept of universality more concrete, 
^such a fixed point is known as non-critical or high temperature 
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Universality classes are characterised by their eigenvalues, A„, (and more generally 
by the dimensions of the operator spectrum and their symmetries) and thus the 
eigenvalues calculated from field theory may, at least in principle, be compared 
with measurements taken in the laboratory. We interpret universality as arising 
due to the fact that near the critical surface the correlation length becomes large. 
As a result the behaviour of the system under consideration is to a large extent 
unaffected by the detailed structure of the model under consideration [9] . Thus 
models 'converge' in the critical domain leading to universality. 

1.9 Critical exponents 

Within the Local Potential Approximation we will use the renormalization group 
to compute experimentally measurable critical exponents, which are related to 
the A„. Specifically, we define for the largest positive A 



To calculate these critical exponents we need to consider the behaviour near the 
fixed point. Suppose we are near a fixed point but not on the critical surface. 
The flow of the action will be dominated by the largest relevant eigenvalue, and 
thus by (1.32) we can write 



1 



(1.34) 



and for the negative A closest to zero 



a;=|A|. 



(1.35) 




(1.36) 
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where is the largest eigenvalue and is the eigenperturbation associated 

with it. We see that /l(0) has scaling dimension — and is associated with a 
coupling of dimension Al0. Now consider the situation at A = Aq. There will be 
an eigenperturbation corresponding to the deviation from criticality. This will be 
the dominant eigenperturbation determining the statistical state of the system. 
Therefore, we associate this operator with 

5Vk^ = (T - T,)^(0) + 0((T - T,f). (1.37) 

As we know the coefficient of has scaling dimension Al, we see that the 

scaling of the temperature difference is also A^. Therefore as the correlation 
length ^ has scaling dimension —1, we must have that 

e ~ (T-T,)-^ (1.38) 
= {T-T,)-r (1.39) 

Corrections to this scaling behaviour are given by the subleading exponents, and 
are dominated by uj. Specifically 

e ~ (T - T^)-" + a^{T - T^f''-^^. (1.40) 

Strictly speaking this is only pertinent to a critical point with one relevant direc- 
tion. With more than one relevant eigenvalue one normally considers a universal 
scaling function of the relevant parameters and then again the least negative 
eigenvalue quantifies the first corrections to scaling. 



^recall that A has been used to trade dimensionful coupUngs for dimensionless ones 
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1.10 /^-functions 



Let us remark that the flow equation in (1.28) can be rewritten in terms of a 
complete set of parameters {a.k.a. the coupling constants). In general (away 
from a continuum limit) this will require an infinite number of couplings. Then 
(1.28) may be rewritten in terms of the /3- functions defined as 

P'W) = f . (1.41) 

By the definition (1.30), the /3-functions vanish at a fixed point where g"^ = gl for 
all i. Later we will perform this reformulation within an approximation. 



1.11 Anomalous dimension 

As described it is convenient to rewrite the flow equations in terms of dimen- 
sionless quantities, denoted here with a subscript dim. We do this by scaling the 
couphngs and fields by appropriate powers of A. Hence, for classical scalar field 
theory we would let 

= 0dimA^(^-2) (1.42) 

such that the kinetic term in the action is overall dimensionless. However, to take 
account of wave function renormalization in the corresponding quantum theory 
we use the anomalous dimension, rj, by letting 

= 0<iimA^(^-2+''). (1.43) 

Compare this with (1.23) to deduce that this is indeed the anomalous dimension 
introduced by (1.11). Likewise we scale the other operators (0^* for integer i — 
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0, 1, ■ ■ ■)0 using a set of parameters which we denote ■ji which may be deduced 
from the renormahzation group (for example 0^ ~ ^^^.^A'^^^+^i). Similarly the 
couplings, g^, scale in an 7, dependent way leaving the action dimensionless. 
Then considering the scaling behaviour of the couplings with A we deduce which 
correspond to relevant operators and irrelevant operators. For example, the mass 
is found to scale as 

m'r^ml^A'-'^ (1.44) 

which is clearly relevant provided 71 < 2. In fact, in this way it may be possible for 
classically relevant operators to become irrelevant and visa-versa when quantum 
corrections are included. From (1.44) it is clear that a fixed point with 71 < 2 
corresponds to choices of couplings and interactions in the Lagrangian such that 
it produces a theory which is massless, regardless of the value of the dimensionless 
massQ {rri'^-^), whereas 71 > 2 corresponds to an infinitely massive theory, and 
7i = 2 is marginal. 



1.12 The Momentum expansion 

In practice, exact flow equations prove very difficult to handle. Here we will be 
interested in using the so called Local Potential Approximation (LPA hereafter) 
as a method of deriving flow equations which can be solved with relative ease. 
We consider the effective Lagrangian as an expansion in powers of momentum 
and within the LPA truncate at the zeroth order, hence discarding all momentum 
dependence. For example, for scalar field theory we could write 

S'" = I d'^x [v{<P,t) + \{d,<pf^ (1.45) 

^here we are primarily interested in a Z2-invariant theory and concentrate on terms which 
maintain this symmetry 

^recall that at a fixed point, the dimensionless mass to^^^ is fixed under variations in A 
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in D dimensions. This approximation has been extensively investigated [10-13] 
and has produced some competitive results, as will be discussed. Clearly as all 
quantum corrections to the kinetic term are neglected, the LPA is equivalent to 
assuming vanishing anomalous dimension. 

By including successive momentum dependent terms it is expected that a se- 
quence of converging improvements over the LPA will be generated [14] allowing 
for non-zero anomalous dimension. For example at first order we may parama- 
terise S^^f by 



where a = 1, • • • , iV for an A^-component theory. This leads to three partial 
differential flow equations for V, K and Z. For one component field theory we 
can simplify the effective action by eliminating K and Z in favour of a single 
quantity, K -\- (f)'^Z [14] . We should note, that in general the expansion (beyond 
LPA) is found to destroy 'reparameterisation symmetry' [15], the requirement 
that physics should not depend on the normalisation of (f) {i.e. (j) Q(f) and 
J — > ^). In some cases this problem can be circumvented by a careful choice of 
cutoff. 



The title 'Renormalization Group' is somewhat unfortunate because it may be 
confused with the perturbative (Callen Symanzik) renormalization group that ap- 
pears in most textbook treatments of quantum field theory. It should be stressed 
that more generally the renormalization group is merely a framework, a set of 
ideas which may also be applied to problems quite unrelated to field theory. Thus, 



I 





1.13 Discussion 
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Kadanoff blocking simply represents the traditional introduction to renormaliza- 
tion from a field theory perspective. However, whatever the motivation, these 
methods have the common feature that they lead to mathematical equations de- 
scribing renormalization group flows in some complicated parameter space. It 
is the study of these flows, and what information they yield about the system, 
which is the essence of renormalization group theory. With quantum fleld the- 
ory, the renormalization group methods become particularly advantageous when 
we consider the continuum limit. We have found that a fixed point of the flow 
equation corresponds to the continuum limit. Thus, by solving for a flxed point 
we are solving directly in the continuum. 

In this chapter and throughout the rest of this thesis we will be primarily inter- 
ested in calculating critical exponents, which provide a convenient testing ground 
for the Local Potential Approximation. However, the ultimate motivation for a 
flxed point search runs much deeper. To make a quantum fleld theory well deflned, 
traditionally we impose an ultra-violet cutoff, A^. Then, the bare parameters in 
the Lagrangian are chosen functions of Aq, such that as the cutoff is removed, 
we obtain a theory with only flnite mass scales. On the natural scale Aq this 
means that all the mass scales must vanish as Aq — > oo. This in turn implies 
that the choice of bare parameters must be such that it results in a fixed point 
(as discussed in section 1.8). The Standard Model of particle physics is based 
entirely around the Gaussian flxed point, which has problems. The scalar flelds 
of the Higgs sector may deflne a trivial theory about the Gaussian flxed point, 
in the sense that as Aq is removed {i.e. Aq oo), all choices of bare couplings 
result in a free theory. Also, there has been no explanation for why the masses 
and couplings take the values observed experimentally. One possibility is that the 
Standard Model might be deflned about a non-trivial flxed point. Such a flxed 
point would have couplings at some non-zero values and would be an interacting 
scale free theory (conformal fleld theory) . 
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Chapter 2 



Scalar field theory 



In the following analysis we make more concrete some of the appealing argu- 
ments and concepts presented in chapter one. As discussed, the ultimate goal is 
to derive and solve flow equations for non-trivial theories. However, it is useful 
to begin with the simplest scalar field theory to develop techniques and under- 
standing. As indicated, we will use critical exponents as a testing ground for 
the approximation scheme under consideration. In the first part of this chapter, 
we investigate standard methods {i.e. truncations and shooting) as applied to 
a Z2-invariant one component field theory. In particular we will be interested 
in the non-trivial Wilson-Fisher fixed point found in three dimensions, however 
some analysis between two and four dimensions is also presented. Subsequently, 
we develop a new variational procedure and compare with the standard methods 
presented. The large N limit of the Local Potential Approximation is solved 
analytically and found to yield the eigenvalue spectrum exactly. In the closing 
section we compare these results with those due to other leading methods. 
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Figure 2.1: A schematic representation of the momentum cutoffs used 

2.1 The Wilsonian effective action 



We begin by making some more precise statements about the Wilsonian effective 
action [16] discussed in section 1.7. The generating functional 



/■Ao 

Jo 



V(f) e 



—Sba 



is replaced with 



Z= [ V(pe 
Jo 



(2.1) 



(2.2) 



(In practice, we ignore a multiplicative A-dependent factor which does not con- 
tribute in any correlation function). Wc arc then free to lower A as represented 
schematically in figure 2.1. 

Consider an A^-component scalar field theory with an ultra-violet (UV hereafter) 
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1 ^ 

Figure 2.2: An ultra-violet cutoff function 

cutoff imposed using the kinetic term. We write the A-dependent action as 

Sren = ^4>a-^uV '^l^a + Sa (2.3) 

where the dots have the usual interpretation as representing an integral over 
position spaceQ and where a runs from 1 to A^, in momentum space Ajjv = 

2 2 

and Cuvij^) is an UV cutoff function. We insist Cuv is a function of ensuring 
that it is dimensionless and Lorentz invariant, and that it falls off sufficiently 
rapidly such that the theory is efficiently regulated. Furthermore, we require 
that Cuv is a profile that acts as a UV cutoff, i.e. C[/y(0) = 1 and Cuv for 
g ^ oo. It might take the form shown in figure 2.2, however it sometimes proves 
convenient to choose Cuv as a smooth [14] {e.g., power law) function. Note that 
as Cuv ^ for g > A in figure 2.2, Sren in (2.3) diverges providing the required 
cutoff for the partition function in (2.2). Similarly we introduce an infra-red (IR 
hereafter) cutoff by modifying the propagator ^ to Ajr = where Cir{j2) is 
a profile with the properties C/r(0) = and Cjr ^ 1 as g ^ oo. We require 

^e.g. f.g = fxQxy = / d^xf{x)g{x,y) and tr{g) ^ g^x ^ J d^xg{x,x) 
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that the two cutoffs are related by 



CiR { 1^) + Cvv (l^ ) = 1- (2.4) 



No source terms have been included in the generating functional, (2.1) and (2.2). 
If a source term, normally denoted J, is included in the bare action, then the 
renormalization procedure will result in additional J dependent terms in the 
renormalised action. In this case the correlation functions, given by 

' ' ' Z[.n (2.5) 



Z[0] 5J{x) 6J{y) 



equal to <0(j;)0(y)> at the bare level will be completely unaltered by the renor- 
malization procedure. 



2.2 The Polchinski flow equation 

The Polchinski flow equation [17] governs the flow of Sa under varying A. Here 
we present a constructive proof following the methods found in ref. [18]. The 
partition function for component scalar field theory with propagator A~^, 
arbitrary bare action [0] and source J, can be written 

where the field (f) has a flavour index a which runs from 1 to A^. Consider the 
following integral which we will show reduces to (2.6). These fields 4>ir and 4>uv 
are integrated over freely. 

Z[J] = j'D(l)jj{D(f)^y e~^'^uv^uv-<t>uv-h't'iR-^7if't>iR-SAoi't'iR+'t>uv1 + M<t>ia+4>uv) (2.7) 
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where 

Aip) = Aj^ip) + A^yip) (2.8) 

and 

4>a = 4>IR + 4>UV- (2.9) 

For simplicity of notation we have suppressed the flavour index a on the right side 
of (2.9). Physically, we interpret the flelds as the momentum modes higher 
than and the (j)^y flelds as the modes that are lower than A. By changing 
variables in (2.7), (f)j^ = 4> — 4>uv and (j)uy = (p^y + A[/yA~^0, and performing 
a Gaussian integral we reproduce the partition function in (2.6) up to a constant 
of proportionality. 

Now consider integrating only over the higher modes: 

ZaI-JAuv] = Jvcpjji e-^''''^'-^^^-'^'^-''^0'f'iR+'^uv]+M<PiH+M, (2.10) 

Using (2.9) to eliminate (pjj^ in favour of (p, it is straightforward to show that Z\ 
does not depend on both J and (|)^y independently in the sense that 

Za[J, (j)uv] = e^-^-'^'R-^+^-'^uv-SAl^iR-J+'f'uv] (2.11) 

for some functional S\ which depends only on the sum A jr. J + (t>uv- If we 
restrict the support of J to low energy modes only, i.e. we set Ajj^.J = 0, (2.11) 
simplifles. However from (2.10) and (2.7), 

Z[J] = jv^^y ZA[J,4>uy]e-"^^UV^uUuV^ (2.12) 

and we see that Sa is nothing more than the interaction part of the Wilsonian 
effective action, as deflned in (2.3). The exact flow equation follows readily from 

^i.e. > A^, 



25 



the fact that (2.10) depends on A only through the (j)jj^.Aj^.(j)jj^ term. Thus, 
differentiating Z/^ with respect to A we immediately obtain the flow equation for 
Za: 

^Z44.av. = ^uv) . (^A,-) . (A _ z,. (2.13) 

Substituting (2.11), performing a change of variables $ = Ajr.J + (^uv and 
relabelling = $, yields Polchinski's version of the Wilson flow equation [17]: 

OSa 15Sa dAuv SSa 1 fdAuv S'Sa\ 



dA 2 (50 ■ 9A ■ 50 2 V ' 

where the trace represents an integral over position space. In terms of the cutoff 
function, C = Cuy, this is written 



OSk 1 / \ SSa ^, 



tr C'.-^ - I^.C:-^ . (2.15) 



dk A3 V V 5<Pa5<i)aJ 50a 50a 

where we have reintroduced the flavour index. 



2.3 The Legendre flow equation 

The Legendre flow equation will prove particularly useful for us as it can be solved 
exactly in the large limit. However, we will extensively use the Legendre flow 
equation for Fermions and thus omit a calculation of the large A^ limit for scalar 
field theory^ Consider the action with kinetic term modified to include an infra- 
red (IR) cutoff, as described in section 2.1. Then the generating functional can 
be written 

Z[J]= yp0 e-5'^-^«-'^''-^A+-^-'^- (2.16) 
■^we will however demonstrate that the LPA provides exact exponents in this limit 
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Figure 2.3: An infra-red cutoff function 

where a runs from 1 to as before, but here A/r = where Cm might take 
the form shown in figure 2.3. We begin by defining the generating functional of 
connected Greens functions by W[J] = ln(Z[J]) and introducing the Legendre 
effective action, rA[0"^], 



+ -0:.A7^.0^ = -W[J] + J„.0^, 



(2.17) 



where in the classical limit /i — > 0, Fa — > -Sa. The classical field is defined as 



<t>'a 



5W 



(2.18) 



and corresponds to the expectation value of the field. Differentiating (2.17) twice 
with respect to 0'^, and then separately twice with respect to J, we derive 



_fw_ 

6JaSJb 



-1 



(2.19) 
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We also differentiate (2.17) once with respect to A, to find 



dA 



1 , 9A7^^, dW 



(2.20) 



Finally we find an expression for the right side of (2.20) by differentiating the 
partition function with respect to A, 

dW ^ _1JW_ dAjl6W_l fdAjl fW_\ 

OA 2 6J' dA '6J 2 \ dA ' 5J5J ) ^ ' ' 

which then leads to the exact Legendre flow equation upon substituting (2.19), 



aA 2 V 

where 

Kb = + Kb^ll (2.23) 

The equations are best appreciated graphically. Wc view (2.21) as depicted in 
figure 2.4, which shows how the n-point functions Wnipi, ■ ■ ■ ,Pn) of Wa evolve. 
The vertices are drawn as open circles, whereas the black dots represent the 
two point function, In the case of a sharp cutoff these represent a delta 

function restricting momenta to g = A. The Polchinski flow equation (2.14) 
can be viewed in a similar manner. In fact the Polchinski and Legendre flow 
equations are intimately linked [18]. In particular the Wilsonian action behaves 
hke the generator of the connected diagrams and these effective actions are thus 
found to be related by a generalised Legendre transform, 

Sa^ = + - 0)a.A7i.(0^ - 0)„. (2.24) 

Eliminating Fa from (2.22) above, we retrieve the Polchinski flow equation (2.14), 
if we impose (2.4). A more detailed discussion of this relationship can be found 
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Figure 2.4: The flow equations for tlie vertices of tlie generating functional Wa_ 
in the literature [18]. 



2.4 Application of the Local Potential Approx- 
imation 

The Local Potential Approximation (LPA) can be apphed with considerable suc- 
cess to cither the Polchinski or Legcndrc flow equation. Here we choose to focus 
our attention on the Polchinski flow equation [17]. Consider the flow equation 
in terms of the cutoff function, C, given by (2.15). Here we impose the LPA by 
writing Sren as (1.45), thus 



.c. 



Jd'^xd^'y V'{cf>{x))C',^V'{cf>{y)) 



(2.25) 



Ha 



Ha 



Fourier transforming. 




(2.26) 



29 



= C'(-^yix-y). (2.27) 

Considering C"(— ^) as an expansion in powers of (— f^) we discard all but the 
leading term, C"(0). Then we find that (2.25) can be written 

^.C'.^ = C'{0) fd^x V''(<l>(x)). (2.28) 

0(Pa 0(pa J 

Later we will assume that C"(0) < 0. In principle the C"(0) could well be taken 
to vanish however we take C"(0) < for the approximation to proceed. This may 
be considered physically unreasonable however the flow equation is ultimately 
independent of C"(0) , suggesting that this assumption may actually be sensible. 
Similarly, we simplify the first term in (2.15) using (2.26) to find 



' 5(t)a{x)5(t)a{y) 



Id'^xd^'yC'—^^^—- (2.29) 
AA'^Jd^'x y"(0(x)) (2.30) 



where we define the dimensionless A through 



^jd%C'(^^)l <0. (2.31) 



(27r) 

Then the 'reduced Polchinski equation' reads 

^ = AA^-V" - ^V^'^. (2.32) 

Scahng to dimensionless quantities via — > 0A^(^~^)-y/[]4[ (where as mentioned 
in section 1.12 77 = in the LPA) and V VK.^ \cf(S))\ using renormalization 
time as defined by (1-29) we finally arrive at a cutoff independent, and hence 
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scheme independent, flow equation, which will be the subject of our attention: 

dV 1 dV fdV\'^ d'^V 

DV--(D-2)4..— -[—]+-,. (2.33) 



dt 2 



a 



The physical predictions (critical exponents) are independent of the scaling^ of 
and V , although the form of the potential will be altered. 

Given the definition of a fixed point a.sV = V^ such that ^ = 0, it is immediately 
noticed that there are two trivial solutions, the Gaussian and High Temperature 
fixed points which are given by 

= (2.34) 



and 



Vr =\<P'-^ (2-35) 



respectively. We will shortly justify the name High Temperature for the latter. 
In addition we will find at least one non-trivial fixed point below four dimensions. 
Prom (2.33) we can deduce that with the exception of the Gaussian solution, in 



the large regime, K(0) ~ 



To compute the critical exponents we perturb about the fixed point solution as 
described in chapter one. We substitute ^(0, t) = K(0) + ^5V{(t),t) into (2.33) 
and work to 0(e). Then by separation of variables we find that the potential 
takes the form of (1.31) as claimed and that the /„(</>) satisfy 

^ = n-D)U + -{D-2Ua^ + 2^^. (2.36) 

For the trivial fixed points we can easily solve to extract the A„. Restricting to 
= 1, we observe that the solutions take the form of a polynomial, fn{<P) = 
0^ + B{A)(f)^^^ + C{A)(f)'^-^ H , with A = 2n for integer n such that /„(0) has 



A<p and V BV for arbitrary constants A and B 
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no divergences. The A„ can then be written in terms of A. In particular, for the 
trivial fixed points we find 

X^ = D- n{D - 2) (2.37) 

and 

X^^ = D-n{D + 2). (2.38) 

Note that the eigenvalues for the Gaussian fixed point are the classical dimension 
of the couplings. We observe, that for the Gaussian fixed point in any dimension, 
the leading exponent defined by (1.34) isQ = |. However, for the High Temper- 
ature fixed point and for all D > we find that all the A„ < implying that there 
are just irrelevant directions around this fixed point. In addition, the exponents 
of the High Temperature fixed point are independent of corresponding to an 
infinitely massive theory with no propagation. 



2.5 Flow and functions 

For illustration, we consider the simple case of the Gaussian fixed point perturbed 
by the mass operator, for a single scalar field. Thus, using the notation introduced 
in section 1.10, we set (for = 1) V{(f),t) = g\t) + ^g^{t)(j)'^. The /3-functions 
for and g^ then follow easily by substitution in (2.33): 

P\g) = ^ = Dg' + g' (2.39) 

and 

/3\g) = ^ = 2g\l-g'). (2.40) 
^this result is in direct agreement with mean field theory [5] 
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Figure 2.5: The /3-function for the mass term of three dimensional one component 
Z2-invariant scalar field theory 

These can be solved for g^{t) and g^{t). As expected, at the Gaussian {g^ = g^ = 
0) and the High Temperature {g^ — and g^ — 1) fixed points, the /9-functions 
vanish. Now we can plot the /3- function for the mass (figure 2.5). We deduce 
that if we add a small mass term to the Gaussian we will flow into the High 
Temperature fixed point. In other words the mass term corresponds to a relevant 
perturbation for the Gaussian fixed point, but an irrelevant perturbation for the 
High Temperature fixed point. 



2.6 Truncations 



We are interested in finding non-trivial fixed points. Here, we outline the method 
of truncations for one component theory. In this approximation method (in 
addition to the LPA) we expand K(0) in powers of (p and truncate the series 
at some finite power, 2p (p = |, 1, |, • • •). Begin by considering as an infinite 
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series in powers of 



(2.41) 



n=0 



such that the fixed point version of (2.33) yields 



- 2)K - DVn = (n + l)(n + 2)K+2 - uvV^V, 



(2.42) 



u+v=n+2 



as the coefficients of a general power, 0". In the last term u > 1 and f > 1 
are integers. We will consider a Z2-invariant theory (even V(0)) which requires 



14 = for odd n. 



Similarly we write an eigenperturbation as a series^, 



n=0 



such that the coefficients of f/)"" in (2.36) yields 



(2.43) 



n 



(D-2)/„-(D-A)/„ = (n + l)(n + 2)/„+2-2 5^ uvVuU (2.44) 



M+D=n+2 



Given that we truncate the series at a power 2p we find that the eigenvalues are 
found by solving the determinant of a (p + 1) x (p + 1) matrix, 



f D-X 2 

2 - 81^2 - A 



det 








-16^4 

-2We 




12 

4 - D - 161^2 - A 
-32^4 



(2.45) 



For illustration consider the simplest case of p = 2, which leads to three simulta- 



^the /„ defined here are distinct to the eigenfunctions used in the latter part of section 2.4 
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neous equations for the V^: 



and 



H = -^V2, (2.46) 



V2 = - 6^4, (2.47) 



(4 - D)V4 = I6V2V4. (2.48) 



Solving, we reproduce both the Gaussian and High Temperature fixed points. 
Additionally we find a non-trivial fixed point, 

4-D (2.49) 

16 ^ 768 ^ ^ ^ 

Immediately we observe that it becomes degenerate with the Gaussian in 4 di- 
mensions and the High Temperature in —4 dimensions. The former is observed 
generally. The first three eigenvalues for the Gaussian and High Temperature 
fixed points, given by (2.37) and (2.38), are reproduced exactly. In principle the 
non-trivial fixed point in three dimensions would have A = 3, | ± | \J^, however 
(2.49) is actually unbounded from below. This is unphysical and will be seen to 
be an artefact of the truncation. 

Generally, the method of truncations proves problematic. The exponents calcu- 
lated in this case are found to initially converge but then fail to converge further 
with increasing p [19]. Most notably spurious solutions are generated for p > 2, 
and some procedure must be sought to distinguish the true fixed points. In some 
cases it has been found that some improvement can be made by expanding around 
the minimum of the potential [20]. The truncation method has been found to 
yield imaginary exponents (a; here) which would imply exotic flows near the fixed 
point. However, later we will find that the exponents are expected to be real, 
excluding the possibility of exotic fiows. 



35 



Figure 2.6: Shooting: (f)c ^ s 

2.7 Shooting 

Exact non-trivial fixed points (within the LPA) are found by solving (2.33) nu- 
merically. We use a method known as shooting [13]. We see that there is a one 
parameter (denoted s) set of solutions, corresponding to a second order differen- 
tial equation with a single boundary condition. For a given value of s we integrate 
out for V{(j)) until it diverges, generally at finite = 0c- Then we 'scan' over 
s and look for non-divergent solutions at s = s^=, which correspond to our fixed 
pointsQ. For consistency, we consider again a Z2-invariant one component theory 
such that V is even in and then set our boundary conditions as 

V{0) = (2.50) 

'^this is a requirement of Griffiths analyticity combined with the required analyticity at the 
origin 
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Figure 2.7: The Wilson-Fisher fixed point 

and 

V'{Q) = (2.51) 

where the s = V"{0) and prime denotes difi'erentiation with respect to (f). The 
first condition is a translation directly from (2.33) and the second is due to the 
symmetry of the theory. A plot of 0c against s is shown for three dimensions, in 
figure 2.6. The fixed points occur for diverging as expected from (2.34) and 
(2.35) the Gaussian occurs at sf = and the High Temperature occurs at sf^ — 
1. The non-trivial Wilson-Fisher fixed point [16] is observed at = —0.229. 
In figure 2.7, we plot V^^. 

Varying the dimension, we observe that there are no non-trivial fixed points at 
or above four dimensions. In exactly four dimensions we find the Wilson-Fisher 
fixed point is degenerate with the Gaussian, where the Gaussian has a marginal 
eigenvalue, A2 = 0. With decreasing dimension, the Wilson- Fisher fixed point 
occurs for decreasing s, such that it moves away from the Gaussian (left in figure 
2.6), and new non-trivial fixed points are created which are initially degener- 
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Figure 2.8: A plot of s* v D for the Wilson-Fisher and two subleading fixed points 

ate with the Gaussian for dimensions where the = (n = 2,3, ■ ■ ■). These 
marginal eigenvalues allow the addition of the corresponding perturbation to the 
Gaussian, creating a new fixed point as previously discussed. The dimensions are 
given by 

D = H(!i±il. (2.52) 
n 

Hence, below four dimensions we see at least one non-trivial fixed point whereas 
below three dimensions we see at least two (Appendix 2 A) . At and below two di- 
mensions we observe a continuous line of 'oscillating' fixed points. The behaviour 
of with dimension D is illustrated for the Wilson-Fisher and two subleading 
fixed points in figure 2.8. 

Finally (2.36) is solved for the exponents in a similar manner. Here we impose 
two conditions, 

/(O) = 1 (2.53) 
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and 

/'(O) = (2.54) 

where the first is an arbitrary normahsation and the second is an imposition of 
our choosing (we are restricting to even perturbations). The solutions generally 
diverge and are found to depend on the eigenvalue, A. Thus, similar to before we 
sweep through a range of A and look for acceptable {i.e. non-divergent) /(0). For 
the Wilson-Fisher fixed point, the largest values yield u — 0.6496 and u; — 0.6557. 
Later, these will be compared to other leading methods. 



2.8 Variation 

We begin by introducing p(0, i) = e"^^^'*) and G(0) = e~2('°-2)</>2 ^^^^ ^^^^ 
Polchinski equation (2.33) can be written [21,22] 

where 

- a^JdU G + (1 -2lnp)\ (2.56) 



and a > is a normalisation factor for the (I) measure to be determined later. 

Sp 



Then the fixed point is defined by ^ = 0. An interesting reparameterisation of 



(2.56) is discussed in Appendix 2B. 

As discussed in chapter one we let g'^{t) be a complete set of parameters (the 
coupling constants) for V. Away from fixed points these will be infinite in number. 
Then using (2.55), we write 

fjN J. 
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dt dg^ 



(2.59) 



where (5^ are the /3-f unctions given by (1.41) and M^j is defined by 



N 




d'^ct) G 



dp dp 



(2.60) 



dg '- dg^ 



Perturbing about the fixed point solution by writing g^ = gi + 5g^ = gl + Ve^^, 
(2.59) yields 



such that (assuming that V is not a null vector) we calculate the eigenvalues, A, 
by solving 



The fixed point equation, ^ = 0, suggests the possibility of approximating p by 
a variational ansatz, p = f{(f), (?') where i runs from 1 to a finite M and / is of our 
choosing. Then interpreting the functional derivative in (2.55) to include only 
variations in this restricted set, we arrive again at (2.59) where however, indices 
only run from 1 to M. Thus geometrically, we restrict the fiows to a submanifold 
M. parametrised by g^, - ■ ■ , g^. 

It is important that these methods reproduce at least the trivial fixed points 
accurately. Here we will rederive the Gaussian and High Temperature fixed points 
with the leading exponent using the methods described above. We take V — 




(2.61) 



dg'^dg^ 



(2.62) 




(2.63) 



40 



g^{t) + ^g'^{t)(f)'^ as an ansatz and minimise the functional in (2.56) 



dg' 



a 



N 



(2.64) 



= 



(2.65) 



to arrive at a pair of simultaneous equations for g^ and g^. Performing the 
resulting Gaussian integrals (Appendix 2C), we arrive at the following conditions 



Solving these leads directly to the Gaussian {g^ = g^ = 0) and High Temperature 
[g^ = —j^ and g"^ = 1) solutions. Then differentiating (2.64) with respect to g^ 
and substituting (2.60) into (2.63) and again performing the resultant Gaussian 
integrals we reproduce the leading exponents for the trivial fixed points. For the 
Gaussian (2.63) leads to {X — D){X — 2) — and for the High Temperature we 
find (A-D)(A + 2) = 0. 

2.9 The Wilson-Fisher fixed point by Variation 

The true test of the variation method lies in the accuracy with which it can repro- 
duce the Wilson-Fisher fixed point [22] . Thus, working again with a Z2-invariant 
one component theory and using the simplest non-trivial ansatz, V^(0, t) — g^{t) + 
\g'^{t)(tP' -\- g^{t)(f)^, we find an approximate solution to the Wilson-Fisher fixed 
point, corresponding to gl = 0.05479, gl = -0.13488 and g^ = 0.00773. In figure 
2.9, we plot the resulting form of p* and compare it to the exact non-trivial fixed 




(2.66) 



and 




(2.67) 
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point solution^, solved by shooting. Solving (2.63), we obtain, apart from A = 3, 
V = 0.6347 from the positive eigenvalue, and u = 0.6093 from the remaining 
eigenvalue. These are 2% and 8% away from the exact^ values respectively. 

Clearly, it may be possible to find 'spurious fixed points' that do not well ap- 
proximate the exact solutions, but these will not be stable under changing or 
improving the ansatz manifold Ai. In fact as Ad improves, such fixed points (if 
indeed there are any) will disappear. These characteristics are in marked contrast 
to the general situation for truncations, as discussed earlier. To improve A4 we 
may simply extend the series to include a (p^ term, or alternatively appeal to 
(2.33) to build a more sophisticated ansatz. In practice such parameterisations 
are increasingly difficult to solve and the results presented here are of sufficient 
accuracy. As applied here, the variational method provides no real advantage 
over shooting. The true potential of this method lies in the relative ease in which 
approximations for global flows may be solved and approximate solutions found 
for more than one invariant. 



2.10 The large N limit of the LPA 

The large limit of the LPA to the Polchinski equation can be solved analytically 
and compared with known exact results. It is thus a good testing ground for the 
LPA. In fact, we will find that the LPA yields exact results in this limit giving us 
added confidence in the reliability of this approximation. Using other methods 
[1] various authors have found that, for large A^, there exists a fixed point in 
dimensions 2 < D < 4 with z/ = ■jjz2- Here we reproduce these results, using the 
methods of Morris et. al. [23]. 

^exact within the LPA 
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Figure 2.9: The simplest polynomial variational approximation to the Wilson- 
Fisher fixed point (dashed line) compared to the exact solution (full line) 

We begin by defining z{x) = (f)l{x) > such that within the LPA V can be 
written as a function of z and t only, corresponding to an 0{N) invariant theory. 
Then the Polchinski equation (2.33) can be written 



We observe = and V = the trivial fixed points, however concentrate 

on non-trivial solutions. By considering figure 2.10 (due to (70^) and requiring a 
non-trivial finite limit, g ~ Ng"^ which implies that g jj. Then ~ -^0^ 
which implies cfp' N and thus V ^ N where the mass is independent. Thus 
changing variables so as to scale out the N dependences, and taking the infinite 
N limit the second order derivative in (2.68) vanishes. Then, differentiating with 




(2.68) 



43 



~g ~Ng 

Figure 2.10: A^-dependence of the two leading contributions to the four point 
function 

respect to z and defining W — ^ we arrive at the following flow equation, 

dW dW 

-— + ((D -2)z + 8zW -2)^ ^2W -4W^. (2.69) 
at oz 

Alternatively, shifting to a stationary point in V , Zo-, (thus excluding the High 
Temperature solution) and defining ^ = + we arrive at [23] 



W - ZoW + {D-2) {zo + y)W'-2W^ 2W' -8{zo + y)WW' - AW^. (2.70) 

Here the dot denotes differentiation with respect to t and the prime denotes 
differentiation with respect to y. 

However, because we have shifted to a minimum we require W{0, t) — Oin (2.70), 
which yields W'{0,t) = and thus W{y,t) = (by expanding in powers of y), 
unless 

Zo ^{D- 2)zo - 2. (2.71) 
Substituting this non-Gaussian solution back into (2.70) yields 

W+{D- 2)yW' -2W = -8{zo + y)WW' - 4W^. (2.72) 
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Then, expanding in powers of y, 

oo 

W{y,t) = J2Wn{t)y\ (2.73) 

n=l 

yields a series of equations, 

Wi + {{D-2)-2)Wi = -8zoW^, (2.74) 
W2 + (2(D - 2) - 2)W2 = -2AzoWiW2-9Wl, (2.75) 
H/„+(n(D-2)-2)W„ = -^{n+l)zoW^W^---. (2.76) 

The final equation holds for n > 1 and the dots stand for terms containing 
products of W^'s with m < n. Thus the Wn{t) are all soluble in terms of Wm, 
with m < n, and Zo{t), and truncations to some finite n are exact. 

Now consider a fixed point, W{y,t) — W{y) and Zo{t) — Zq. Then by (2.71) 

however the original field squared it must be positive, requiring D > 2. 

Note that since there is only one solution for Zq, W{y) can only cross the axis 
once. Substituting this into (2.74) yields Wi — Q (the Gaussian) and 

Since W{y) crosses the y-axis only once, we must have W\ > Q which implies 
we must have D < A (as we have established D > 2), so that W{y) > for 
y > 0. Otherwise the potential is unbounded below. Substituting (2.78) into the 
series of equations (2.74) to (2.76) we see that the W„ exist and are unique {e.g. 
Wo = —^W?) 

2 8-5D 1/ 
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We deduce the leading exponent by perturbing the Zo in (2.71) about the fixed 
point solution, Zo = ^o* + 8zo, yielding 5zo = {D — 2)6zo, such that there is an 
eigenperturbation Szo oc e'^°*, with Aq = -D — 2 implying 

jy = 2 < D < 4. (2.79) 

Now we can systematically generate the subleading exponents by perturbing 
about fixed point solutions in (2.74) to (2.76). For example, by writing Wi — 
Wu + SWi in (2.74) (with 5zo = 0), we arrive at 5Wi = (D - A)5Wi, or 
SWi oc e(^-^)*, yielding \i = D - 4. Then in (2.75) we find a new perturba- 
tion which takes the form 5zo{t) — 0, 5Wi{t) — and 5W2{t) oc e(^~^)*, yielding 
A2 = -D — 6. In this way we generate the set of eigenvalues in the large N hmit, 
An = D-2nforn = l,2, 



2.11 Discussion 



Referring to the numerical methods, the most outstanding feature is that they 
allow us to exhaustively search the entire infinite dimensional space using a sin- 
gle parameter and to make physical predictions for complex systems which can 
be measured experimentally. The results for u and lu found by shooting and 
the variational method compare favourably with other leading methods (table 
2.1) [1]. It is not yet known how to apply the momentum expansion (beyond 
LPA) to the Polchinski fiow equation in a reparameterisation invariant [15] way. 
However the momentum expansion at O(p^) has been applied to the Legendre 
fiow equation and shows a modest improvement in the exponents, over the LPA 
[14] . Traditionally it is difficult to estimate the errors introduced using these ap- 
proximation schemes, other than by comparison with the other methods, though 
recently progress has been made [14]. The LPA is found to be competitive, at 
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least for small N and away from D — 2 where rj becomes non-negligible. 

The success of large N limit of the LPA to the Polchinski equation is somewhat 
fortunate. It should be noted that the flow equation (2.68) is not itself exact, 
but belongs to a large class of equations which yield the correct exponents [23]. 
Nevertheless, it is an important result and this limit will be the focus of our 
attention for Fermionic field theory. It is hoped that solving this limit will lead 
to an increased understanding of the problem and thus help us to deal with finite 
N theories. Similarly, truncations have been found to be useful, albeit in a limited 
capacity. In fact from (2.62) we can now see that provided the potential V is 
real, we expect the eigenvalues, A„, to be also real, as claimed. This contrasts 
with the imaginary exponents which occur using the method of truncations. 



Method 


V 


UJ 


Lattice calculation 


0.6305 




e-expansion at O(e^) 


0.6310 


0.81 


Six loop perturbation series 


0.6300 


0.79 


Local potential approximation (Pol) 


0.6496 


0.6557 


LPA Variation method (Pol) 


0.6347 


0.6093 


Local potential approximation (Leg) 


0.6604 


0.6285 


Momentum expansion at 0{p'^) (Leg) 


0.620 


0.898 



Table 2.1: Exponents for three dimensional one component Z2-invariant scalar 
field theory 

The above comparisons are made using results found in [1] and [14]. 



47 



Appendices 



Appendix 2A: Non-trivial fixed points 

We include the data for non-trivial fixed points below four dimensions, including 
the first two subleading fixed points (table 2.2). 



D 










Wilson-Fisher 


1st subleading 


2nd subleading 


4.00 


o.ooot 






3.90 


-0.009 






3.80 


-0.021 






3.70 


-0.035 






3.60 


-0.051 






3.50 


-0.070 






3.40 


-0.092 






3.30 


-0.119 






3.20 


-0.149 






3.10 


-0.186 






3.00 


-0.229 


o.ooot 
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2.90 


-0.280 


0.005 


- 


2.80 


-0.342 


0.015 


- 


2.70 


-0.417 


0.032 


- 


2.60 


-0.512 


0.062 


-0.001 


2.50 


-0.632 


0.108 


-0.008 


2.45 


-0.707 


0.139 


-0.015 


2.40 


-0.793 


0.178 


-0.028 


2.35 


-0.895 


0.226 


-0.049 


2.30 


-1.019 


0.284 


-0.082 


2.25 


-1.172 


0.355 


-0.137 


2.20 


-1.368 


0.440 


-0.226 


2.15 


-1.633 


0.542 


-0.378 


2.10 


-2.027 


0.661 


-0.651 


2.05 


-2.738 


0.811 


-1.249 


2.00 


-oo 


1.000^ 


-oo 



Table 2.2: Wilson-Fisher and subleading fixed points. 



t Fixed points degenerate with the Gaussian. 

I Fixed points degenerate with the High Temperature. 

Appendix 2B: Functional reparameterisation 

We can reparametcrisc the functional, in such a manner, that the first term 
in (2.56) looks like a kinetic term. As a consequence the second term can be 
interpreted as a time varying potential, unless we deal with two dimensions. 
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when it is static. This is done by making the following change of variables, 

ci0 = (ir e3(2-^)<^'. (2.80) 

Then, replacing T with —S we arrive at the following, with the desired kinetic 
term and time varying potential: 

S = ldr(^^p'-f{T)p'(^\np-^-^^ (2.81) 

where 

fir) = J e^(^-^)^^ (2.82) 
and the dot denotes differentiation with respect to r. 

Now we explicitly view (2.81) as an action, and think of a particle at position p, 
at time|^ r. From (2.82) it is easy to see that in dimension greater than two, /(t), 
and thus the potential, decreases with time, whereas for dimension less than two 
the reverse is true. The form of the potential, 

V(p) = p2(^lnp-0 (2.83) 

is shown in figure 2.11. 

We now view fixed points as motions in the potential well which are either at 
rest permanently, or, only come to rest at p = in an infinite time (recall that 
in the large regime, K(0) ~ |0^)- example, the Gaussian fixed point 
corresponds to a particle which simply remains static in the minimum of the 
well, whereas the High Temperature fixed point corresponds to a solution which 
lies at p = exp(l) ^ 1.40 at r = and then moves towards p = as r ^ oo. 
The solution for the non-trivial fixed point has been found (figure 2.9), and this 

^not to be confused with the renormalization time defined by (1.29) 
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Figure 2.11: Time varying potential in reparameterised functional 

corresponds to a particle which hes at p 0.92 at r = 0. For this solution the 
particle oscillates once before coming to rest at p = 0, again as r ^ oo. These 
motions can be deduced directly by observing the behaviour of p with cf). 



Appendix 2C: Gaussian integration 



For completeness, we list the set of Gaussian integrals used: 



J — ( 



dx e 2"^^ 




27r 
a 



(2.84) 



found by considering 7q and changing variables, (x, y) (r, 6). Then we define 



-f 



(2.85) 
(2.86) 



6=0 
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and calculate 



= -J-, (2.87) 

a\ a 

h = (2-88) 
V a 

= ^y^. (2.89) 
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Chapter 3 



Zamalodchikov's C-function 



In the preceding chapter we argued that the Local Potential Approximation pro- 
vides competitive estimates of the critical exponents, at least for scalar field the- 
ory. The theoretical analysis presented here further vindicates and increases our 
confidence in the approximation. We begin by introducing conformal field theory, 
leaving a more detailed analysis to reference [4] . We review Zamalodchikov's C- 
theorem [24] using the introductory material. The C-theorem, as presented here, 
only makes sense in two dimensions, although a number of groups have sought 
to generalise these ideas to other dimensions [25-36]. However, within the Local 
Potential Approximation we construct a function, C, which has an appropriate 
generalisation of the Zamalodchikov properties, in general dimension, D. Our 
C-function 'counts' degrees of freedom at fixed points. We normalise such that 
it counts one per Gaussian scalar and zero at the High Temperature fixed point 
(corresponding to an infinitely massive theory). 
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3.1 The confer mal group 



As discussed in chapter one, we are interested in non-trivial fixed points which 
correspond to an interacting scale free theory, known as a conformal field the- 
ory. Here we will work mostly in two dimensions, however recent work [22,25] 
has extended some of the aspects presented to four dimensions. We begin by 
considering a D-dimcnsional Ricmannian manifold upon which a metric, Q/^i,, is 
defined. The interval between two points on the manifold is given by 

ds'^ = g^.^dx^dx" (3.1) 

and a coordinate transformation, x — > x', generates a change in the metric tensor 
given by 

9i^Ax) ^ 9'i.u{x') = Q^-Q^9ap{x). (3.2) 

The conformal group is defined to be the subgroup of the transformations that 
leaves the metric unchanged apart from a non-zero differential scaling function 
fl{x) such that 

9,xu{x) g'^^{x') = n{x)gf,^{x). (3.3) 

These metrics arc said to be conformally related; the angle between vectors and 
the ratios of magnitudes of vectors are the same on both manifolds. However, 
since the length of the vectors differ, this transformation corresponds to a rescal- 
ing. The Poincare group leaves the metric totally invariant and is thus a subgroup 
of the conformal group with Q = 1. 
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3.2 Special properties in two dimensions 



The infinitesimal generators of tlie group can be determined by tlie analysis of 
an infinitesimal coordinate transformation, — > x'^ = — e^{x), which causes 
a line element change. In fiat space ((7^^ = 

ds'^ ^ ds"^ = dx ^dx^ (3.4) 
= ds^ - (d^e^ix) + d^e^{x))dx^dx'', (3.5) 

where we have used (3.2) and defined 9^ = Note that if e'^(x) = e'^, i.e. a 
translation, then the line element is invariant under such a transformation. For 
the transformation to be conformal we require by (3.3) and (3.5) that 

n{x)ds^ = ds^ - (9^e^ + d,,ey)dx^dxr (3.6) 

Then using (3.1) we arrive at 

5^^{l - n{x))dx''dx'' = {d^e^ + d^e^)dx''dx'' (3.7) 

and for general dx^ 

5^,{l-^l{x)) = d^e, + d,e^. (3.8) 
Then, using S^^S^p = S'^ = D, we find that the scaling function is given by 

n{x) = l-^id.e). (3.9) 

where d.e = d'^e^. 

In two dimensions, substituting (3.9) into (3.8) yields the Cauchy-Riemann con- 
^ working in Euclidean space 
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ditions for the analj^icity of sq + iei , 

doto = diti and doei = —dieo. (3.10) 

This motivates us to use complex coordinates defined by 

z — + ix^ and z = x'^ — ix^, (3-11) 

and write the infinitesimal transformation as 

€{z) = eo + iei and e{z) — eo — iei. (3-12) 

Then the conformal transformations coincide with the coordinate transformations 

z^f{z) and z ^ f{z) (3.13) 

which have an infinite dimensional local algebra, as f{z) and f{z) can be any 
analytic functions. These are the generators of the local conformal algebra with 
Euchdean line element ds^ — {dx^Y + {dx^Y — dzdz transforming as 

2 

dzdz = fldzdz. (3-14) 
However, if we take the transformed coordinates 

z' ^{1- enZ^)z (3.15) 

and 

z^z' ^{l- enZ^)z (3.16) 
where n ranges over the integers, then it is easily verified that the infinitesimal 



ds' 



dl 
dz 
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generators are given by 



Ln^-z^+'—^-z^+% (3.17) 

and 

K = (3.18) 
Now we can construct the algebra: 

[Lm: K] = {m- n)L^+n, (3.19) 

[Lm, -^n] = (m - n)Lm+n (3.20) 

and 

[L„,L„]=0. (3.21) 

The latter condition, (3.21), shows that the local conformal algebra is a direct 
sum of two independent subalgebras. The action of the D — 2 conformal group 
thus factorises into independent actions on z and z and the Greens functions of 
the theory are such that z and z can be treated as independent variables. 



3.3 The stress tensor as the generator of scahng 

Let us now consider a Euclidean action S{(f)) depending only on the field and its 
first partial derivatives, invariant under translation (associated with momentum 
conservation), rotation (angular momentum conservation) and dilation. Let us 
perform an infinitesimal change of variables: 

x'a^ Xfj.- ^u.{x). (3.22) 
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As discussed in Appendix 3A, due to translational invariance, the variation of 
the action involves only the partial derivatives of e^(a;): 

55(0) = Jd''xT^,{x)d^e,{x) (3.23) 

where T^,/(a;) is the stress tensor. Imposing other conformal symmetries, we 
deduce some required properties of the stress tensor [1]. For example rotation 
invariance implies 5S vanishes for 

= Kux"" (3.24) 

in which A^j, is an arbitrary antisymmetric matrix. Therefore the stress tensor 
must be symmetric: 

Tfjii, — Tufj,- (3.25) 
Similarly, dilation invariance corresponds to 

= ^^tJ. (3-26) 
and implies the vanishing of the trace of the stress tensor: 

T^^ = e = 0. (3.27) 

In complex coordinates we have d^z — d^z — Q as z and z are treated independent. 
We refer the metric to the complex coordinate frame using (3.2), 



I ° J I , (3.28) 
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and similarly calculate the stress tensor, 

T.. = ^(Too-2zTio-Tn), (3.29) 

T,, = ^{Too + 2iTw - Tu) (3.30) 

and 

Tzz — = '^T^^. (3.31) 

Furthermore, we can calculate the derivatives of (3.29) and (3.30) with respect to 
z and z respectively, such that by conservation of the stress tensor (see Appendix 
3A), d^T^y = 0, we deduce 

dzTzz + dzT,, = dzT + ^dzQ = (3.32) 

and similarly 

dzT-z-z + dzTz-z = dzT + Us = 0. (3.33) 

Here we have defined T{z,z) = Tzz{z,z) and T{z,z) = Tzz{z,z), 9 is given by 
(3.27) and we have used (3.31). Finally we note that for a scale invariant field 
theory^ (3.31) vanishes by (3.27) and we find d^T = dzT = which implies 
that the stress tensor factorises into independent pieces, T{z, z) = T{z) and 
f{z,z) = fiz). 



^corresponding to a fixed point 
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3.4 Conformal weights 



In general, the fields in a conformal theory obey the two dimensional equivalent 
of the tensor transformation law, 

where {h,h) are the conformal weights of We observe that T = T^z and T = T^z 
have conformal weights (2,0) and (0,2) respectively. Similarly the trace of the 
stress tensor, © = ^Tzz, has conformal weight (1,1). Now consider two fields 
^i{z,z) and $2(-2,^) with conformal weights (/ii,/ii) and (/i2,^2) respectively. 
Then by (3.34) < $i(z, z)$2(0, 0) >^ aC^i+MaCMa) < $i(/(^), /(^))$2(0, 0) > 
for /(2;) = az. For a rotation a = e*^, and thus we require 

<t,(.,^)t.(0,0)>=-;JgL_ (3.35) 

for the transformation to have no effect. In (3.35) E is an arbitrary scalar func- 
tion. If, in addition we have scale invariance (corresponding to a fixed point), i.e. 
a — re*^, we have a stronger condition for invariance, 

E 

< ^Az, ^)$2(0, 0) >= , , r r (3.36) 
where now E must be a constant. 
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z 



z 




z 



w 



Figure 3.1: Evaluation of the commutator (3.38) on the conformal plane 

3.5 The operator product expansion 

To generate the operator product expansion we construct a conserved charge on 
the conformal 2;-plane, 



associated with the conserved stress tensor discussed above. For the full confor- 
mal group there are an infinite set of conserved currents (associated with scale 
invariance at a fixed point), e{z)T{z) and their antiholomorphic partners. These 
charges can then be used to generate an infinitesimal symmetry variation in the 
field, $, given by 



where we have combined the two z integrations in the commutator into a sin- 
gle contour integral around w [4] (figure 3.1) and normalised with a factor 
However, taking^ f{z) = z + t{z) and using (3.34) we can expand for small e to 

■^consistent with the transformation used in section 3.2 for smah e(z) = —z"'^^ 




(3.37) 




(3.38) 



(3.39) 
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find 



<^{w,w) {l + hd,e){l + hd^e){<^> + ed,<^> + eds<^) (3.40) 
= $ + h{d,e)^ + h{d,e)^ + ed,^ + e9,-$. (3.41) 

Then, for (3.39) to be consistent with (3.41) we deduce that we must have 

T(z)$(w, w) = - — ^-^$(w, w) + -^—^^^w, w) + --- (3.42) 

[Z — W)'' z — w 

and similarly 

f{z)^{w, w) = -—^^<^iw, w) + -^d^^{w, w) + --- (3.43) 

[z — wy z — w 

where the dots contain contributions analytic in z. We may verify (3.42) and 
(3.43) by substituting into (3.39) and performing the contour integration, repro- 
ducing the variation in (3.41). The expressions (3.42) and (3.43) are known as 
the short distance operator expansion and apply generally to a tensor of weight 
{h,h). 

Alternatively, in the corresponding quantum theory, T{z) is written as a normal 
ordered product of two fields and Wicks rules can be used to reproduce the 
operator expansions above. Then, it is found that in the case of the expansion 
of T with itself an additional term arises due to double Wick contractions, which 
can be written [4] 

T{z)T{w) = ^ + —^T{w) + -^^^d^T{w) + ■■■ (3.44) 
2[z — w)* [z — wY [z — wY 

where c is known as the conformal anomaly^ [37], and the coefficient of T{w) 
reflects the fact that T has a conformal weight h = 2. A similar expression for 

^also known as the Virasoro central charge 
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the expansion of T with itself can be found, however we will now concentrate 
purely on the unbarred case. Observe that as the stress tensor written as a 
normal ordered product annihilates the vacuum, the two-point function is given 

by 

< T{z)T{Q) >= (3.45) 

and thus it is expected that c > for a unitary theory. Note that this is precisely 
the form required at a fixed point by (3.36). 



3.6 The Virasoro algebra 

To generate the Virasoro algebra [4] it is useful to define a Laurent expansion of 
the stress tensor in terms of operator modes, 

T{z) = Y,z-^-^K (3.46) 

neZ 

with 

K^-^.jdzz^+^T{z). (3.47) 

Classically these are simply Fourier modes but become operators modes in the 
context of quantum field theory. From (3.15) we substitute e[z) — —z^~^^ into 
(3.39) to verify that (3.47) is indeed the infinitesimal generator. The exponent, 
z~'^~^, is chosen to give the operators, L„, a scaling dimension of n under a scaling 
z az, consistent with (3.18), since then T[z) — ^ ^T{az) by (3.44) and thus 
Ln ^ oI'Lr, in (3.47). 

To derive the algebra satisfying the above modes we must commute the operators 
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Ln to give a double contour integral, 

li.. - f / ^ / S - / S / ^1 ."«T(.)u,"'«r(.). (3.48) 



27ri / 27ri J 27ri / 27ri / 

The z contour integration can then be performed by initially fixing w and deform- 
ing the difference between the two z integrations in (3.48) into a single contour 
around w (figure 3.2). If we then use the stress tensor self product, (3.44), and 
finally integrate over w we can write the commutator as. 



2{z-w)'^ 



+ ^^-^T(«;) + ^^^a^T(«;) + • • • ) (3.49) 

[z — wY [z — wy ' 



,.n-2 



= ([ —w"'^' ( — {n + l)n(n-l)w' 
J 2m Vl2 

+ 2{n + l)w"^(^/;) + w^'+^d^.Tiw) + •••). (3.50) 



where we have explicitly performed the z integration. Re-expressing the third 
integral in (3.50) by parts and combining it with the second term, we finally 
arrive at 

[Ln, Lm] = ^{n'-n)j 1^^"+"^"' + {n - m) j ^w^+^^'T{w). (3.51) 

Referring to (3.47) we identify the second integral as Lm+n and then performing 
the first integral we derive the Virasoro algebra. 



Lm] = —{n^ - n)5n+m,o + (n - m)Ln+m (3.52) 



and similarly 



[Ln, Lm] = - n)6n+m,0 + {u - m)Ln+m- (3.53) 
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• 0' 

Figure 3.2: Evaluation of the commutator (3.48) on the conformal plane 

Every conformally invariant quantum field theory has a particular representation 
of this algebra with some value of c or c. We observe that the Virasoro algebra 
reduces to the classical case forc = c = 0, (3.19) and (3.20). 

3.7 Counting 

In the following sections of this chapter we motivate, review, and discuss the 
implications of Zamalodchikov's celebrated C-theorem [24] for two dimensional 
quantum field theory, and establish three important properties for his C-function. 
The third property depends on the conformal anomaly and thus only makes sense 
in two dimensions. A number of groups have sought to generalise these ideas to 
higher dimensions [25]. Then, within the LPA, we will display a C-function 
which has the first two of these properties and an appropriate generalisation 
of the third property in any dimension D. At fixed points, our C-function is 
extensive, videlicet additive in mutually non-interacting degrees of freedom, as is 
also true of the conformal anomaly. 

Suppose that the fields of the theory form two non-interacting sets. Let us write 

0a = 0a^'' when the field belongs to the first set, and (pa = (ptf^ when it belongs to 
the second set. Similarly, as they are non-interacting, the couplings g split into 
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two sets denoted g^^^ and g^"^^ and the action can be written 

S[(f),t] = + (3.54) 

Then, the variation in the action and thus the stress tensor, can be written 
as a sum of two separate parts as can be seen from (3.23). Furthermore, the 
generators given by (3.47) separate into two independent pieces and thus the 
Virasoro algebra in (3.52) leads to the conclusion that c is extensive at fixed 
points, i.e. 

c{g.)=c{gi'^)+cigi'^). (3.55) 
This counting will be the third property of our C-function derived below. 

The interpretation given to this is one of counting degrees of freedom. At the 
Gaussian fixed point, = 0, the conformal anomaly counts one degree of 
freedom per scalar. At the High Temperature fixed point, V^^ = — ^, the 
conformal anomaly vanishes, corresponding to an infinitely massive theory]^ with 
no degrees of freedom. Thus, we will normalise such that our C counts one for 
each Gaussian scalar and zero at the High Temperature fixed point. 

3.8 Zamalodchikov's C-theorem 

Zamalodchikov's C-thereom provides an explicit geometric framework for the 
space of two dimensional quantum field theories [21,22] and demonstrates the 
irreversibility of renormalization group flows, thus proving that exotic flows are 
missing. Zamalodchikov established three properties for his C-function, 

^in units of A 
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I) There exists a function C{g) > of such a nature that 

where form an infinite set of parameters (the couphng constants discussed in 
section 1.10) and the beta functions are defined by (1-41). 

II) C{g) is stationary only at the fixed points of the renormalization group, i.e. 
when g{t) = g^. 

III) The value of C{g) at the fixed point g^ is the same as the corresponding 
conformal anomaly [37] . (This property thus only makes sense in two dimensions.) 

The proof relies on rotational invariance, reflection positivity and the conservation 
of the stress tensor (a property which is a general consequence of translational 
invariance and is also valid away from the critical point). Consider some point 
on a renormalization group trajectory (refer to figure 1.4) specified by a set of 
couplings, g\ For the time being we suppress the dependence on g^. As discussed, 
away from the fixed point the Hamiltonian is no longer invariant under dilations 
and in addition to the components T and T the stress tensor has non-zero trace, 
©. These components have conformal weights (2,0), (0,2) and (1,1) respectively, 
so that by (3.35) we can write their two-point functions as 

<Tiz,z)TiO,0)> = (3.57) 
<e{z,z)T{0,0)> = <T{z,z)e{0,0)> = (3.58) 

Z'^Z 

H(zz) 

< 6(^,^)6(0,0) > = (3.59) 

z z 

where F, G and H are non-trivial scalar functions. Taking the correlation function 
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of (3.32) with T(0, 0) and 0(0, 0) yields two equations 




(3.60) 



(3.61) 



where we define F — zzF'{zz). Then, defining 



C ^2F -G --H 
8 



(3.62) 



and ehminating G from (3.60) and (3.61) yields 



(3.63) 



Refiection positivity requires that < 0© > > 0, so that by (3.59) H > 0. Then 



Now imagine making a renormalization group transformation such that the lattice 
spacing a a + 5a. Since our C-f unction is dimensionless, this is equivalent to 



If we now define C{g) = C{ro,g^), where Tq is some arbitrary but fixed length 
scale, then given that ^ < 0, we deduce the first Zamalodchikov property for 
C. Moreover, C is stationary if and only ii H = 0, which by refiection positivity, 
implies 6 = 0, so that the theory is scale invariant (condition for (3.27)) and 
thus corresponds to a fixed point, which implies property two. Finally at such a 
fixed point, the vanishing of the trace of the stress tensor (©) by (3.27) implies 



from (3.63) C < which imphes C" < as ^2; > 0. This in turn implies that C 
is a non-increasing function of r = (zz)^ — ((a;°)^ + (x^)^)^ > as C" = 2r^. 



sending r ^ r — 6r, and the coupling constants will flow according to the 
renormalization group equation. 




(3.64) 
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that G = // = by (3.58) and (3.59). Thus (3.62) becomes F = |C, so that 
comparing (3.45) with (3.57) shows that C — c verifying the final Zamalodchikov 
property. 

3.9 Mij as a metric 

We begin by noticing that My, defined by (2.60), can be interpreted as a positive 
definite metric [22] in the space of quantum field theories. Initially we consider 
the interval, 



The final observation (positive definiteness) holds providing V is real (corre- 
sponding to a unitary theory in Minkowski space). The equality in (3.67) is only 
reached if dg'^ = 0. Perhaps more importantly we observe that the interval is 
indeed invariant under a change of coordinates. 



(3.65) 




(3.66) 



> 0. 



(3.67) 



(3.68) 




dg^dg^ 



(3.69) 



= Mkidg^dgK 



(3.70) 
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3.10 A C-function representation of the LPA 



We define our C(5i)-function through 



(3.71) 



where b > 1 and is defined by (2.56). This C-function satisfies our appropriate 
generahsation of Zamalodchikov's properties in general dimension, D. Generally 

^^^*-lnW. (3.72) 

dg' dg' ^ ' ^ ' 

However, by differentiating (3.71) with respect to 51* and using (2.59) we find 

dh'^ 4 

and therefore 

^ = (3-74) 
where we have defined Mij through M^j — {J^nb)Mij. Thus, we find that 

dt dg' ^ ^ 

= -/3'Mij/3^. (3.76) 

Using an argument similar to that used to derive (3.67) we can now write 

< (3.78) 

confirming the first Zamalodchikov property. Given that the /9-functions vanish 
at a fixed point, the second property follows immediately from (3.74). 
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The appropriate generalisation of the third property is found by supposing that 
the fields form two mutually non-interacting sets. Using our previous notation 

we write (pa = 4>a'^ when the field belongs to the first set, and (pa = (p^a^ when 
it belongs to the second set, such that the couplings g also split into two sets 
denoted g'^^^ and g'^'^^ and by (3.54) the potential can be written 

y ((/), t) = ((/)(^) , t) + 1/(2) (0(2) , t) (3.79) 

and thus p factorises: p = p^^^p^'^\ However, substituting the fixed point equation 
^ = into (2.56), we derive that fixed points p{(j),t) = p*(0) satisfy 

nP*]-ja''jd''cl>Gpl (3.80) 

and we immediately notice that J^[p*\ — ^J^[p2'^]J^[p^^^]. Thus from (3.71) we 
observe that at fixed points our C is extensive, as required by (3.55). 

In two dimensions, at the Gaussian fixed point, the conformal anomaly counts 
one degree of freedom per scalar, i. e. is equal to A^, and at the High Temperature 
fixed point, it vanishes. We will take this to be true in general dimension D. Thus 
the normalisation factors a and b can be uniquely determined by requiring that 
our C agree with this counting in any dimension D. Substituting C = in (3.71), 
and p* = exp(— y/^^) in (3.80) and performing the Gaussian integral yields 

Similarly, for the Gaussian fixed point we substitute C = in (3.71) and = 
exp{-Vf) = 1 in (3.80) and find 

& = eV^y|±|. (3.82) 
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Note that as required, 6 > 1, at least for D > 2. 



3.11 Examples 

Consider the simple example of the Gaussian fixed point perturbed by the mass 
operator, for a sing le scalar field. Thus we set V{(j),t) = g^{t) + lg^{t)(j)'^. Then 
the /5-functions are given by (2.39) and (2.40). Integrating (2.40) yields 



and then (2.39) yields 



g\t) = gy^ + e-^J^^^dr (3.84) 

1 rl .,D/2-l 





where u — e^*"^"^ and k and gl — g^{Q) are constants of integration. Considering 
the limit t — oo, we see that this solution indeed emanates from the Gaussian 
fixed point. We normalised the special solution in (3.85) so that with gl = 0, it 
tends to V^^ as t ^ +oo. 

Since p is Gaussian, is readily determined using (2.56), 
In particular, using (3.81) we verify 
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and 



T'''^ = ^ (3.88) 
4 



for the Gaussian [g^ — g"^ — 0) and High Temperature {g^ — and g'^ — 1) 
respectively. Then using (3.71) and (3.82) we deduce that C"^ — 1 and C^^ — 
as expected. Indeed we may further verify that when = 0, C{t) flows from 1 
to as t runs from — oo to +cxo. Note that if g^ ^ 0, then C{t) ±oo as t ^ oo, 
depending on the sign of g^. This seems in contradiction with the idea that C 
counts degrees of freedom, since evidently the vacuum energy should not figure 
per se in this counting. However, we recall from (3.55) that C is extensive {i.e. 
counts) only at fixed points, and with Qq ^ the system never reaches another 
fixed point as t ^ oo. 

From the Gaussian form of G, wc recognise that at both the Gaussian and High 
Temperature fixed points, the metric (2.60) is diagonalised by choosing the oper- 
ators Oi — ^ to be products of Hermite polynomials in the 0^. Since these 
also turn out to diagonalise g^iggj (g*), they are the eigenperturbations, and the 
corresponding eigenvalues follow straightforwardly. Choosing the Gaussian fixed 
point p* = 1 for example, and again specialising to the case of a single scalar field 
for simplicity, we thus take On — Hn{^^/D — 2), n = 0, 1, • • •. From (2.56) we 
then obtain 



d'^T f dOdO- 



dg^dg^ 

However, using H'^ — 2ni?„_i where prime denotes differentiation with respect 
to the argument, we deduce = n\/D — 20n-i so that (3.89) yields 

. = -DMij + ij{D - 2)Mi_i (3.90) 



dg^dg- 

Then, given that the metric has non-zero components [38] M„„ = 2"+^n!ay^^^, 
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from (2.62) we recover the Gaussian spectrum of eigenvalues, 



A. 



1 d'^T 



(3.91) 




(3.92) 



as expected. The High Temperature fixed point can be solved in a similar manner. 



The C-theorem has the interpretation that renormalization group flows go 'down- 
hill'. As described it rules out the existence of exotic flows and in particular re- 
stricts the possible fixed points to which unstable directions at a given fixed point 
may flow. Various authors [39] have reported on the appealing interpretation of 
the C-function as a kind of entropy of information about the critical system. 
Under renormalization, information is lost about the short distance behaviour of 
the correlation functions, corresponding to a decreasing C . 

The value of the C-function at a fixed point is interpreted as the number of degrees 
of freedom, which may initially seem counter intuitive given that the C-function 
we have defined (3.71) is non-integer at the Wilson-Fisher fixed pointQ. However, 
the following example may convince us that this is in fact not as unusual as we 
might imagine. Consider a ball attached to a solid wall on a spring. Initially 
the ball might have one degree of freedom (position in one dimension). Then 
imagine increasing the spring stiffness in a smooth manner until eventually the 
ball is permanently stationary. The degrees of freedom of the ball have reduced 
smoothly from one to zero. 

^the exact LPA solution yields c = 0.9867 whereas variation yields c = 0.9896 



3.12 Discussion 
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It is probably not possible within the Local Potential Approximation to establish 
a more concrete link between our C-function and Zamalodchikov's C. Of course, 
it would be very interesting to understand if these observations generalise to 
higher orders in the derivative expansion [40] (which likely would allow a direct 
comparison with Zamalodchikov's C), or indeed generalise to an exact expression 
along the present lines. 
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Appendices 



Appendix 3A: The stress tensor 

Generally, if the Lagrangian density L depends only on the field and its 
derivative 9^0(a;), the variation in the Lagrangian can be written 

where we have used 

bd^<^{x) = d^b4>{x). (3.95) 

To derive the equations of motion we assume that the arbitrary variations, b(\){x\ 
vanish at infinity so that the last term in (3.94) integrates to zero in the variation 
of the action. However we suppose that 0(x) satisfies the equations of motion 
and take (50(x) not as an arbitrary variation, but a translation: 

b<^{x) = -e''d^4){x). (3.96) 
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Then the first bracket in (3.94) vanishes and substituting (3.96) we find 



However, £ is not invariant under a global transformation, rather we have 

5jC = -e^'di.C. (3.98) 

Thus, from (3.97) and (3.98) 

which, since it is true for all e'^, implies 

d^T^, = (3.100) 

where we have defined 

T^.^j^M-Si:C. (3.101) 

The stress tensor, T^,^, is the conserved quantity corresponding to translational 
invariance. 

Now consider the variation of the action with a local translation e'^(x), 

SS = Jd^xSC (3.102) 

where we have used (3.95) and (3.96). Replacing the inner brackets of the first 

77 



term in (3.104) with di,C and integrating the last term by parts we finally arrive 
at 

5S = Jd^xe''{x)d,(^j^^^d,(l){x)-6i^C^ (3.105) 
= Jd^xe''{x)d^'Ti,^ (3.106) 

by the definition of the stress tensor (3.101). 
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Chapter 4 



Fermionic field theory 



In the previous two chapters we have assembled a considerable portfolio of ev- 
idence in support of the Local Potential Approximation. However, it must be 
admitted that (here and in general) most success has arisen from focusing on 
purely Bosonic field theories, and ultimately we must extend these constructions 
to include Fermions in a reliable manner. Here we will be specifically interested 
in the Fermionic Polchinski and Legendre flow equations. For a finite number of 
components, the LPA to the Legendre flow equation is significantly more difficult 
to implement than for the scalars and the Polchinski flow equation will form the 
focus of our attention. However, a technical problem involving chiral symme- 
try emerges and we find that the Polchinski LPA is trivial. These results are 
compared with exact solutions obtained in the large N limit using the Fermionic 
Legendre flow equation^. As a result we develop an understanding of the Polchin- 
ski LPA for a finite number of components, which are presented towards the end 
of this chapter. 



only the Legendre flow equation is exact in the large N limit [23] 
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4.1 Introduction 



In this chapter we will be interested in Fermionic field theory, and in particular 
spinors, which we label ipaa, where here a runs over internal spinor indices from 
1 to Ni and a is the flavour index running from 1 to Nf. These satisfy the well 
known anti-commutation relations including 

{^aa, i^pb} = and {ijjaa, ^/'^'^J = h6/6abS{x - y). (4.1) 

This so called Grassmann nature implies that V-'iV-'i = where i runs over all 
indices, and thus that a function of a finite number of Fermion fields can be 
expanded as an exact finite series. It is useful to use the traceless gamma ma- 
trices [1] {ix — 0, •••,3) and ^^^^^^^^^^s^ which satisfy {7^,7^^} = and 
(7°)^ = J where / is the unit matrix. The adjoint spinor is then defined by 
^« = V^t/37°^". We will be using derivatives -~ and ^ which both act on the left 
in the conventional manner. 

A theory is said to be chirally symmetric if it is unchanged under 

V' ^ 75V' (4.2) 

and 

^ ^ V575- (4.3) 

In particular we observe that a chirally invariant theory has no mass term, or 
indeed any odd power of ■^V'. However the kinetic term is acceptable, where 
we define ^ = 7^^^. Ultimately we arc interested in the Gross-Neveu model [41] 
which can be defined by the bare Lagrangian 

C = C(^)^a%a(^) + 9{ra{x)^aa{x)f (4.4) 
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which is seen to be chirally invariant. 



4.2 Conditions on Fermionic potentials 

To derive the condition that must be satisfied by the Fermionic potential we 
impose that the Minkowski Hamiltonian density, 

H[^P{x),,p{x)] = ^{x)M^) - CMmx),i,ix)] (4.5) 

where the Lagrangian is given by 

CMm^):^l^{x)] = iC(^)^a%a(^) - V{^|;\x)^oi^{x)) (4.6) 

and 

n°(a;) = ■ ^ = -#(a;), (4.7) 

dib (x) \ \ / 

must be positive definite. For convenience we suppress the indices on the argu- 
ments of functions. The dots in (4.5) and (4.7) denote differentiation with respect 
to t = By substituting (4.6) and (4.7) into (4.5) we find an alternative ex- 
pression for the Hamiltonian density 

H[i^{x)Mx)] = Vl"(2^) (a.p)/ V(^)/3a + V{i^\x)^^i:{x)) (4.8) 
where a* = 7^7* and = —id^. 

We begin by substituting this Hamiltonian into the equation of motion 



,(x), jd''ynmy),^(y)] 



(4.9) 
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such that using (4.1) we can perform the integrals over delta functions, yielding 



(4.10) 



Comparing this with (4.8) we observe that 



7i[V'(x),^(x)] = #I"(x)V'a„(x) - Vl"(a;, 



+ y(V'^(x)7V(x)). 



(4.11) 



Now consider a term in the potential which takes the form (^1^(a;)7°'?/^(a;))". For 
this term the commutator in (4.11) can be expanded as 



i^aaix), d^y{i^\y)j'i^{y)r 



n d''y[iP^a{x),i^\y)j'i^{y)] 



X (V'ny)7V(y))"-' + o(n^). 



(4.12) 



Using (4.1) we compute the commutator on the right side of (4.12) to find 



i^aaix), d^y{i^\yWi^{y)r 



n^T{x)^'f,PpSx) 



X (V'^(x)7V(x))"-^ + o(n') 



(4.13) 



and thus deduce that the Hamiltonian (4.11) can be written 

n[ij{x)Mx)] = t^i^{x)ij^,{x) - ^r(x)7°.%„(a;)V^'(^n^)7V(^)) 

+ y(V'^(x)7V(a^)) + Oifi'). (4.14) 

Defining z{x) — '0a (^)'0aa(^) ~ V'l"(^)7°a^'0(^)/3a ^® ii^^s require that the quan- 
tity 

V{z)-zV'{z) (4.15) 
must be bounded from below in mean field situations, where hF' can be neglected. 
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More generally we see from (4.12) that it corresponds to neglecting terms involv- 
ing higher derivatives of the potential, i.e. those of the form 




(4.16) 



z 



dz"- 



with n > 1. We will find that these terms may be neglected in the large N 
limit. Furthermore, we will learn that the correct interpretation (at large N) is 
to treat the z{x) as a real variable. At finite N we find that treating z{x) as even 
Grassmann (and thus a potential, V{z), as a finite series) yields nonsense results 
typical of truncations [19,20]. Thus, throughout we treat z{x) as a real number 
and not even Grassmann as one might naively expect. Finally, the condition 
that (4.15) is bounded from below implies that any potential growing larger than 
linear for large z{x), must be bounded from below. 

4.3 The Legendre flow equation 

The structure of Fermionic field theory leads to a Legendre flow equation that 
looks somewhat different to the analogous equation for the scalars (2.22). How- 
ever, in common with scalar theory it can be solved exactly and analytically in the 
large limit. The application of the LPA to Legendre flow equation is non-trivial 
and thus at flnite N we concentrate on the Polchinski formalism. Ultimately we 
will use the large N limit of the Legendre equation to develop an understanding 
of the LPA to Polchinski's equation. Begin by consider the generating functional 
regulated by an infra-red cutoff, similar to scalar field theory. 




A 



-1 /3 
■IRa 



(4.17) 
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where C and ^ are the sources. We begin by deriving the necessary prehminaries 
to finding the fiow equation. Define Z{(, () — e^^^'^^ and introduce the Legendre 
efi^ective action, Fa, by 



so that taking the classical limit, /i — > 0, 5'a — > Fa. Then, if we differentiate 
(4.18) with respect to ■^'^ and ■0'^ separately, we find 



5T 



A 



and 



A7i/ = - C (4.19) 



aa 



^ + ^-lL'r,a - Caa (4-20) 



where the classical fields are defined by 
and 



5W 

- -Y^ (4-21) 



5W 

similar to the scalar case (2.18). Then, by differentiating (4.19) and (4.20) with 
respect to ■0'^ or ■^"^ and differentiating (4.21) and (4.22) once with respect to C, 
or C we can rewrite the following matrix equation 





<5f 


Sip" 


Sip" 


Sib': 


SC 

Silfi 



I 



Sip^ \ \ IL } \0/ 



(4.23) 



as 



13 



5C5C ^C^C I I SiP"SiP" 

_ _ ] \ S'^Va I A -1 

SCSC SC.sC. / „b c« V ^ 
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/ 
/ 



(4.24) 



ad a 



where / is the unit matrix 

With the prehminaries complete we begin the derivation by differentiating the 
generating functional, (4.17), with respect to the cutoff. A, 



dZ 



S 



dA 6Ca 



OA 5Cf 



(4.25) 



or alternatively, substituting Z = , 



-1 /3 



dA 5Caa dA 5C 



dA SCaa 5Ca 
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(4.26) 



Then we differentiate (4.18) with respect to A and use (4.21) and (4.22) to find 



dW 



dA 



dTA 
dA 



(4.27) 



Equating (4.26) and (4.27) and again using (4.21) and (4.22) we finally arrive at 
our flow equation 



dVK 
dA 



1 





dA 



dA- 



dA 



Sisc 



P J 



(4.28) 



which can be rewritten using (4.24), 



dA 



1 



dA- 



dA 



dA 
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X 



S^Ta S^Ta I A -1 



S'^Ta , A -1 5^Ta 



(4.29) 



which will be our starting point for solving the large N limit exactly. Compare 
and contrast (4.29) with (2.22). 



4.4 The exact large limit 

As discussed, the Legendre flow equation is exactly solvable in the large N limit 
and thus we will use this limit to establish some exact results with which we will 
later compare the LPA (as applied to the Polchinski formalism). We begin by 
assuming that we can write Fa = rA[2;(a;)] for some A [23], where we have defined 
the scalar z{x) = ip2{^)'4^aa{^)- Then differentiating Fa with respect to the fields, 
we observe that the diagonal terms of (4.29) become zero and the off diagonal 
terms take the form 

5a.A7i/ + 5^,5/-^ + C^,.^- (4.30) 

Thus the inverse in (4.29) is of the form (aJ+6)^^ = a~^I—a~^ba~^+a~^ba^^ba~^ — 
■ ■ ■ where a is a flavour singlet, / is the unit matrix in flavour space and b is an 
arbitrary matrix in flavour space. With the trace over flavours in (4.29) this 
imphes that the third term in (4.30) contributes neghgibly for large Nf — Saa. 
Taking this limit, (4.29) becomes 



= ^I + ^7n ■ (4-31) 



OA ' y OA a \ Sz 

If for some value of A, Fa is a functional of z only (as assumed) then + Aj^ 
is also a functional of z only, and thus, given that (4.31) is a flrst order differential 
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equation in A, we have that this form is preserved by the flow, i.e. Fa = ^a[z] 
for all A^. Furthermore, the fact that the flow equation (4.31) now contains 
only zero or one point functions allows for a considerable simplification [23]: 
Set z{x) = z a constant. Then, since the effective action may be written as a 
derivative expansion (to all orders), we see that only a potential term of the form 

Ta= [ d^xV{z,A) (4.32) 



survives. Thus, substituting (4.32) into (4.31) and transforming into momentum 
space 

m - ''''' [ ''-"-'^>™^' {• - ) "7 

where we have used = — A7^^^pA7^ and VL{2t[)^ is the solid angle of a 

{D — l)-sphere. Now we introduce an infra-red cutoff function, Cj/j(i|-), such 
that 

A7ii/ = (« + ^^)"VC^/H/, (4.34) 

where a is an unspecified constant^] and ^jjia ~ ~'^i~^a large momentum. 
Note that this cutoff generally breaks the chiral symmetry (later we will choose 
to restrict the cutoff function and the purpose of the seemingly redundant a will 
become clear). Here, as with the scalars, Cir 1 as g — oo such that the 
physics is independent of the infra-red cutoff at scales much larger than A and 
the theory is well regulated. In terms of the cutoff function (4.33) can now be 
written 

rfV /"°° / 1 



^this has resulted in a restriction to a submanifold in the couphng constant space [23] 
■^not to be confused with the indices 
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X [l+ia + tir'C,^-^] ]. (4.35) 



After scaling into dimensionless quantities, V — > A^V, z — > ^z, g — > Ag and 
a — > Aa, and using renormalization time defined by (1-29) this becomes, 



— = DV + {l-D)z—+tNfQj^ dqq^-'[{at' + i)-VC',^,''Cil' 
X {a + t^)/(^I+{a + ti)-'Cj^^^ y (4.36) 

The cutoff function takes the general form, 

C,«/ = 9{qyf + ih{q')iJ (4.37) 



and 



with 



and 



C/"!^ = b{Q')Sj + ^<q')i/ (4.38) 



9 



(4.39) 



where h,c,g and /i are arbitrary functions of the now dimensionless (f. After 
some manipulation (4.36) can be written 



2 I ^2 a 



X ((2g^/i' + /i)-2i^'^)/(6 + ic^)^^ 
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where prime denotes differentiation with respect to the argument, and 



1 + 



ag — hq^ dV 

Q,2 _|_ q2 



(4.42) 



and 



.ah — g dV 



(4.43) 



+ q^ dz 

We choose specifically h = and g = 6^{q, 1), where 0^{q, 1) is a smooth cutoff 
function, satisfying 9e{q, 1) ~ for g < 1 — e, 9^{q,l) ~ 1 for g > 1 + e and 
< < 1. Ultimately we take the sharp cutoff hmit, 9e{q, 1) — > 0{q — 1) as 
e — > 0, and the flow equation will simplify enormously. Note that the cutoff 
breaks the chiral symmetry unless a — 0. Using = g^tr(7'^) = and writing 
N = NfNi = Nf6^", we arrive at 



(4.44) 



Then using (4.42) and (4.43) and the lemma presented in Appendix 4A [18] we 
take the e — > limit. 



dV dV 



In 



t=o 



(4.45) 



so that we flnally reach the desired flow equation by a change of variables. 



DV ^DV - Nfl [ln{t^)] 



t=0 



(4.46) 



yielding 



dV dV 
-^DVHl-D)z- 



q;2 + 1 



In 1 + 



(4.47) 



In the last step wc scaled out the and VL dependencies via V VLNV and 
z VLNz (note that these scalings justify the neglect of terms of the form (4.16) 
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in (4.15)). Differentiating (4.47) with respect to z and defining 1^ = ^ we find 



The flow equation in its forms (4.47) and (4.48) will form the focus of our at- 
tention. We realise the Gaussian solution, W — V = but in the following will 
concentrate on non-trivial solutions. 

4.5 Solutions in the large N limit 

We now solve the large N limit [42,43] of the Legendre flow equation ((4.47) 
and (4.48)) for the fixed points and eigenvalue spectrum. Initially we choose 
q; = to analyse specifically the solutions with a cutoff introduced in a chirally 
invariant manner, however ultimately we will allow a to be non-zero to make 
meaningful comparisons with the Polchinski approach to be developed later. The 
eigenvalues are found by perturbing the potential in (4.47). We substitute V — 
K + use the fixed point version of (4.47) and expand the logarithmic 

term, yielding to 0(e) 




(4.48) 




(4.49) 



Using the fixed point version of (4.48) it is straightforward to see that 




(A, 



n 



(4.50) 
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and it is easily verified that 



C/n = i^iW^f-^". (4.51) 

where Ki is an arbitrary constant. Insisting that the eigenf unctions, Un, are 
analytic and providing that crosses the z-axis we thus deduce that the power 
in (4.51) must be a positive integer, n, and therefore that 

Xn^D-n (4.52) 

regardless of a. (Note, that these eigenvalues are degenerate with the Gaussian 
in two dimensions.) 

Specialising to three dimensions the exact flow equation, in its form (4.48), can 
be solved analytically using the method of characteristics. The solution is found 
to be 



1 / . 2 / W 



z — —5 — - I a{a — 3)ln 



+ 2(3Q;'-l)tan-^(W^ + Q;)^ ^""'^ 



(a2 + 1)2 

- ^ + /i(^e-)iy^ (4.53) 

CK^ + 1 

where the arbitrary function, /i(H^e~*), is set constant (denoted by /i) for a fixed 
point solution. However, we find we need to impose conditions on the range of 
/i for the solution to be considered physical. We require that the solution W 
is defined for all z and that (4.15), and thus the negative potential, is bounded 
from below. After a little analysis of (4.53) we deduce the large W behaviour as 

a"' + 1 
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Figure 4.1: Solution to the large N limit of the Legendre flow equation in three 
dimensions with cutoff introduced in a chirally symmetric manner 



which leads to a restriction in the range of fi. 

Q!^ + 1 + 1 

Here we have imposed that W{z) > for large negative z and that W{z) < 
for large positive z. Furthermore wc demand that W is not double valued in z. 
This is done by simply differentiating (4.53) and studying its turning points. By 
experiment we find that (4.55) provides a good approximation to this restriction. 

Choosing a — yields a fiow equation for the cutoff introduced in such a manner 
as to maintain chiral symmetry which has the following solution: 

z ^ -2W - 2Whan-\W) + fiW\ (4.56) 

A chirally symmetric potential must be even in z and thus we observe that (4.56) 
actually breaks this symmetry unless we set fi — 0. This solution is plotted in 
figure 4.1. It is noticed that independent of /i, the solution will always satisfy 
W{0) = and thus the eigenvalue spectrum is indeed given by (4.52), An = 3 — n. 
These results are in direct agreement with those of Zinn- Justin [43]. The marginal 
eigenvalue, A3 — 0, is associated with a line of fixed points parameterised by /i. 
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f, 




\ disallowed 



a 



-71 



Figure 4.2: Region of {a, fi) space allowed by the boundedness condition (4.55) 

We find that the double valued condition discussed above restricts the acceptable 
solutions to the range -3.0708 < /i < 3.0708. 

Allowing a to be non-zero, we can study the symmetry breaking scenario. For 
example we could start by fixing fi and altering a, such that for the two bound- 
edness conditions, (4.55), to be consistent we require < i. Thus the region 
of 'parameter space' allowed is constrained (figure 4.2). However, generally the 
double valued condition restricts the region allowed even further. For example 
for a = 0.25 (4.55) implies -2.4024 < /i < 2.4024 whereas in actual fact we find 
—2.2489 < /i < 2.1325. We recognise that (4.53) is unchanged if we simulta- 
neously let a — >■ —a, z —z and W —W (this is discrete chiral symmetry) 
and thus focus on just a > 0. The imposition of (4.55) and the observation that 
(4.53) has no divergences has the implication that all the acceptable solutions 
satisfy W{z) — for some z and thus for all these fixed points the eigenvalue 
spectrum is again A„ = 3 — n. 
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4.6 The Polchinski flow equation 

Here we present a constructive proof of the Polchinski flow equation for a purely 
Fermionic theory. This will form the focus of our attention when we apply the 
LPA. Begin by considering the generating functional regulated by an overall mo- 
mentum cutoff Aq, made explicit later. For a Fermion field ip{x) with propagator 

A(p) and arbitrary bare action S\^['4j,iIj], we write 

Z(C.C) ^ /(#)(#).- ^ * - ^ - (4.67) 

where, as before C, and C, are sources, a and (5 run over internal spinor indices 
from 1 to Nj^ and a is the flavour index which runs from 1 to Nf. Following the 
methods of [18] presented in chapter two we split the propagator and field into 
high and low momentum parts, 

A(p) = A,^(p) + A^^(p) (4.58) 

and 

V' = V'/ii + i^uv (4.59) 
respectively. Then we can write (4.57) as 

Z(C, C) = / [di^uv) (#c7y)e- ^^^"^ "^^^'^ ZSuv^ ^uv, C, C] (4-60) 

where 



(4.61) 
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Figure 4.3: Ultra-violet and Infra-red cutoff functions 

We verify (4.60) by eliminating ■0^^ and 'ipjj^ using (4.59) and substituting ipjjy — 
il^jjy + Ajjy ^jJ to retrieve (4.57) up to a multiplicative factor which we ignore. 
Wc now make the momentum cutoffs explicit by defining Am = [^^(p, A) — 
9eip, Ao)] A and Auv = [1 ~ de{p, A)] A where similar to our previous analysis 9^ is 
a smooth function satisfying 9^{p, A) for p < A — e, ^^(p, A) 1 for p > A — e 
and < ^e(p, A) < 1 (figure 4.3). The derivation of the flow equation starts by 
substituting ipm — — ipuv into (4.61), yielding 



X 



y"(#)(#)e 



= g '^UVa ^IRa '^UV Ha "^Aot^'^l 



(4.62) 



Performing the Gaussian integral in (4.62) 



= e~ ^^^^^ ^^^/^^ ~ ^A„[^4l + ^^uv AJfl+ClS A,^/ [A7^ ^c'^ +^1/3'' (4.63) 



= e 



(4.64) 



so that (4.60) becomes 



Z{U) = Jidi^uvWuv) 



X e ^t^^" + ^" 



'^t/V'a ^Uvf '^UVffa + ''a ^IRcf ^/3a 



(4.65) 



for some functional, S\. 



95 



Now consider that C and C only couple to the low momentum modes, i.e. such 
that (AjR = Ajji( = 0. Then, in this limit (4.65) becomes 

Z{C,C) = hd%y){d%y)e- "^-Va ^U'v! ^UV,. 



X e '^UVa C«a + CS' -^UVoca " ^ f,[i^uV ^'^Uv'\ (4.66) 

and thus S\ coincides with the Wilsonian effective action. Differentiating (4.61) 
with respect to the cutoff. A, gives 

^^+^„,Vi!^7^_,„,) z,. (4.67) 

/ Pa 



OA \SC ^"Va a \5C 



Defining tj^^ = ^c/vaa + ^r/ Cpa and tjj^ = i^uva + Cf A^^^" and substituting 
(4.64) into (4.67) we arrive at our flow equation. 



-eta — a > 



where 



It is possible to confirm this solution by eliminating S\ in favour of Sr 



ijj'^ A^y^ i/j^^ + Sa and then showing that |f = as required. Note the 
similarity between (2.14) and (4.68). 



4.7 Application of the Local Potential Approx- 
imation 

Previous attempts to solve the flow equation [44] (4.68) have involved introducing 

2 

an ultra-violet cutoff function Cuv{j^), as described for the scalars. However, 
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typically various other approximations have been made in addition to the LP7\Q, 
leading to results of dubious reliability. Here we only require that the cutoff 
function be dimensionless and Lorentz invariant, and so for the Fermions we 
allow C = Cuv to be a function of ^, C(z|-), and thus write 

Ac/y/ = -^C^^-\^ (4.70) 

where C(0) = 1 and C falls to zero sufficiently quickly as g ^ oo, such that the 
theory is well regulated. 

Substituting (4.70) into (4.68) and replacing ip with ■?/), 

where the prime denotes differentiation with respect to its argument. Now we 
introduce the LPA by writing 

Sa = J rf^xV^(V',V';A) (4.72) 

such that V does not contain any derivative terms. Similarly we restrict C(^). 
Within this approximation we find that (4.71) reduces to 

where It = A~^{2tt)^^ J d^qC'{^/A) for a dimensionless constant r and the 
C"(0) has been transformed into configuration space. By scaling V — VA^, 
ijj — ^ tljA^^^~^\ ijj 'ijjA^^^~^\ defining the renormalization time by (1.29) 
and finally scaling out the remaining cutoff function dependencies we derive the 



*most notably truncations 
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following scheme independent flow equation: 

dV 1,^ dV ^„dV\ dV dV d'^V 

where all quantities are now dimensionlcss. Note the striking similarity between 
this and the scalar case, (2.33). The Gaussian and High Temperature solutions 
are given by 

K"" = (4.75) 



and 



respectively. 



By perturbing about a fixed point solution, V{ip, ip, t) = V^{ip, ip)+e fni'^P, '0)e'^"*, 
where K is the fixed point under consideration and e is a small parameter, we 
easily find an equation for the perturbations around a fixed point (ignoring 0(e^)), 

Xnfn ^DU-hD-l)( + C 



2^ ' V "-9^- diP'^ 

dlPaadilj^ dlPaadilj^ di/jaaOi/jf ^ ' ^ 

This can then be used to find the exponents which are related to the A„ in the 
familiar way (section 2.4). As for the scalars, the eigenvalues for the Gaussian 
fixed point are given by the mass dimension of the couplings: 

X^^D-n{D- 1) (4.78) 

where n is a positive integer. The eigenvalues for the High Temperature fixed 
point are independent of A'^, 

X^^ = D-n{D + l), (4.79) 



98 



as expected for an infinitely massive theory (seen by calculating the (5 function 
for the mass, similar to that described for scalar field theory). 



4.8 The large N limit of the LPA 

We now consider the large N limit of the LPA Polchinski fiow equation (4.74) 
and compare it to the exact results presented in section 4.5. Similar to before we 
define z{x) = ip2{x)ip^^{x), and consider V — V{z,t). Then the fiow equation 
becomes 

_ = W-(p-l). + 7V)^-.(^) (4.80) 

Scaling out the N dependence via V NV and z Nz and taking the limit, 
the double derivative on the right of (4.80) vanishes. We retrieve the trivial fixed 
points — and V^'^ — ^ + z. Then differentiating (4.80) with respect to z, 

dW o dW 

-W^ + ((l-D)z-2zW -1)^ (4.81) 
at oz 

where W = ^. 

As before we expand the potential in (4.80) as V = K + e Xln UnG^"'^ for small e. 
The usual manipulation leads to an equation for the eigenfunctions 



Un dz (1- D)z-2zW^-l' 

Substituting the trivial solutions into (4.82) we rapidly reproduce the eigenvalue 
spectrum for the Gaussian and High Temperature fixed points. We rewrite (4.82) 

as 

^ - ...... (4.83) 



dW. ^{{l-D)z-2zW,-l) 
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Figure 4.4: Solution to the large N limit of the LPA Polchinski flow equation in 
three dimensions 



(4.84) 



where we have used the fixed point version of (4.81). Integrating (4.84) we solve 
for Un in terms of W^^ 



D-X„ 



(4.85) 



where K2 is a constant of integration. Thus we deduce that if W^, = 1 for some 
z then for [/„ to be analytic at that z, the exponent in (4.85) must be a positive 
integer (n), implying Xn — D — n, whereas if — for some z we would have 
Xn^ D + n. 

Again, solving using the method of characteristics we find the solution in three 
dimensions to be 



1 



(ly- 1)4 \2 



(4.86) 



where /2(^zt) is an arbitrary function. For a fixed point solution /2( ^'Li ) must 
be set to a constant (/2). The f2 — solution is shown in figure 4.4. However, 
from (4.86) we see that for large z, z —-^ and furthermore that we have 
exactly one divergence at 1^ = 1, unless /2 = — |. Therefore by continuity, we 
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Figure 4.5: Non- divergent solution to the large N limit of the LPA Polchinski 
flow equation in three dimensions 

deduce that none of these fixed points can represent acceptable fixed points^. 
The remaining possibility, /2 = — | is plotted in figure 4.5 and is discarded for 
similar reasons. Hence there are no acceptable non-trivial fixed points in three 
dimensions. 



4.9 Finite N 

To study finite theories we use (4.80). We start by substituting a series ansatz, 
truncated at some finite power of z. Generally we find solutions degenerate with 
the Gaussian at D = 2. However, similar to the case with scalar field theory we 
quickly generate large numbers of spurious solutions. The simplest ansatz yields 
the following solution. 



2D \2 A J V8 32 J 

which also collapses to the High Temperature fixed point m D = —2, in a manner 
reminiscent of scalar field theory. However, as discussed we treat z{x) as a real 



^The reader can convince himself or herself of this by graphical analysis 
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number and not as an even Grassmann and therefore we solve (4.80) numerically. 

However, solving numerically, we find that non-trivial fixed points entirely dis- 
appear, unlike with scalar field theory. As discussed in section 2.7 to solve for 
a fixed point we simply set ^ = and integrate (4.80) out from a given V{0) 
until we encounter a divergence at z = Zc. We then invoke Griffiths analyticityQ 
and look for the finite number of non-divergent solutions [13], i.e. ones for which 
Zc —>■ oo. However, due to the fact that the cutoff function is not chirally invari- 
ant, V cannot be assumed even in z. Thus we integrate out for both positive and 
negative z. Figure 4.6 shows the results for three dimensional four component 
theory. For z > and ^(0) < we ffnd that the solution to (4.80) quickly 
diverges whereas all solutions with V{0) > are non-divergent and seemingly 
acceptable ffxed points. However, as discussed, we have to integrate out (4.80) 
for negative z. In this case the solutions diverge for all V{0) with the exception 
of the Gaussian and High Temperature ffxed points indicated. These qualitative 
features do not change with A^. Thus we conclude that consistent with the large 
N limit of the LPA to Polchinski's equation we ffnd no acceptable non-trivial 
ffxed points in three dimensions. New divergences in Zc are generated for z < 
at the Gaussian critical dimensions, those for which the A^, given by (4.78), are 
marginal: 



where n = 1, 2, ■ ■ •. Thus, for D < 2 we ffnd a new divergence in Zc as can be 
seen from ffgure 4.7. This divergence is degenerate with the Gaussian in exactly 
two dimensions, but moves away from the Gaussian as we lower the dimension. 
It is straightforward to see that this does not correspond to a ffxed point as this 
solution diverges in the positive z direction. However for D < | the divergence 
occurs for V{0) > and corresponds to an acceptable ffxed point, labelled NTFP 



^and analyticity at the origin 
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Figure 4.6: Zc versus V{0) for the Fermionic Polchinski flow equation at N 
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Figure 4.7: versus V"(0) for the Fermionic Polchinski flow equation at N 



and -D = I 
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Figure 4.8: Zc versus 1^(0) for the Fermionic Polchinski flow equation at 
and -D = I 
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in figure 4.8. This fixed point becomes degenerate with the High Temperature at 
D = 1, analogous to the first subleading fixed point for scalar field theory. 



4.10 Discussion 

For the Polchinski flow equation, we are at liberty to introduce the cutoff in some 

2 

chirally symmetric manner, for example by insisting that Cuv is a function of 
[44]. However, if we then multiply C^y by the kinetic term in the conventional 
manner, we quickly realise that the right side of the exact flow equation (4.68) 
always carries a ^, and thus, is trivial within the Local Potential Approximation. 
This problem has been alleviated here by allowing Cuv = Cuvij^)^ although it 
should be noted that the non-trivial LPA to Polchinski's equation can be derived, 

2 

albeit less aesthetically, whilst keeping Cuv a function of j2. However, in this 
case, the inverse cutoff function must be introduced in some non-conventional 
manner, such as multiplying 2^, which is then summed with the kinetic termQ. 
Hence we might introduce the cutoff as 



which is still afflicted with a chirally non-invariant term. Thus, in order to obtain 
a non-trivial LPA version of (4.68) we are essentially forced to break chiral sym- 
metry with the cutoff function. To date, there is no known method of including 
momentum dependent terms within an expansion scheme for the Polchinski flow 
equation, without breaking reparameterisation invariance [15]. 

Exact analysis in the large N limit shows that only a small range of symmetry 

^such that the right hand side of the Polchinski equation contains derivative independent 
terms 




(4.89) 
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breaking cutoffs reproduce tfie non-trivial fixed points with the eigenvalue spec- 
trum associated with the symmetry preserving cutoff. Thus it is not surprising 
that non-trivial fixed points are not seen in the LPA to the Polchinski flow equa- 
tion. Generally, we interpret the results obtained from (4.80) and (4.81) as being 
analogous to letting a. become large in (4.47) and (4.48) such that the chirally 
invariant component in the latter case is not dominant. Prom figure 4.2 the line 
of fixed points associated with the large N hmit is seen to vanish for o? > \. 
With this interpretation the results are seen to be entirely consistent. Thus, it 
has to be concluded that the hope inspired by the successes of scalar field theory 
discussed in chapters two and three are not yet fulfilled within Fermionic field 
theories. 
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Appendices 



Appendix 4A: Lemma 

We present the lemma [18] used to (4.45): 

6J{9{p,A),A)^6{A-p) fdtf{t,p) (4.90) 

as e ^ 0, where f{6e,A) is any function whose dependence on the second argu- 
ment (A) remains continuous at A = p in the limit e = 0. This follows from the 
identity 

(4.91) 



5,m(p,A),A)= ||- / dtf(t,A')] 



A'=A 



where here 5e(p,A) — —-^6^(p, A). This can be seen by noting that as e — > 
the lower limit on the integral becomes 1 for p > A and the integral vanishes, 
whereas it becomes for p < A. Thus the integral becomes a step function in A 
but with height dtf{t, A'). Thus differentiating with respect to A reproduces 
(4.90). 
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Chapter 5 



Concluding remarks 



5.1 The Local Potential Approximation 

In the context of scalar field theory, there is little doubt that the Local Poten- 
tial Approximation yields reliable results, comparable with other leading meth- 
ods (table 2.1). However, we have found that the additional approximation of 
truncations, whereby the potential is expanded as a series in the field, 0, and 
terminated at some arbitrary finite power is unreliable. In contrast, standard 
techniques which involve no further approximation, such as shooting, are found 
to be excellent methods of generating the physical quantities of interest, e.g. crit- 
ical exponents. The new variational methods presented here provide an equally 
reliable alternative, with the advantage that we expect approximate solutions of 
LPA's for more than one invariant to be found with relative ease. Furthermore, 
the C-function defined using the LPA to Polchinski's fiow equation provides ad- 
ditional justification and support for this approximation scheme. 

In the case of Fermionic field theory the Legendre fiow equation is significantly 
more difficult to deal with than that for scalar field theory, at least at finite A^. 
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However, within the LPA, the Polchinski flow equation is very similar to the 
corresponding scalar case. Thus, it is unfortunate that the LPA to the Polchinski 
equation fails to produce non-trivial fixed points for Fermionic field theory, due 
to our cutoff function breaking chiral symmetry. It should be considered whether 
re-introduction of scalars might lead to non-trivial fixed points once more. The 
generahsation of (2.33) and (4.74) follows immediately. 



where i is the flavour index for the scalar fields and the notation used for the 
Fermion fields is consistent with chapter four. However the work presented here 
indicates that ultimately it might be necessary to focus on other flow equations, 
such as the Legendre equation. 

5.2 Remaining issues 

Given that the LPA has been found to be competitive, both in terms of simplicity 
of implementation and the numerical estimates of the exponents, at least for scalar 
field theories, it is hoped that ultimately Fermionic field theories will pose nothing 
more than a problem of technique. An interesting question which we have not 
answered is whether there exists some method of extracting results relevant to 
chirally invariant theories from those presented here. We have recently initiated 
some interesting work analogous to that found in ref. [45] which is due to be 
published alongside the work presented in chapter four. 

However, one of the major remaining obstacles is to deal with gauge theories. 
It is found that introducing a momentum cutoff destroys gauge symmetry and 



dV 
'dt 




2 



(5.1) 



a 
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ultimately leads to flow equations with gauge dependent parameters. Clearly, it 
is difficult to make reliable physical predictions when ultimately they are gauge 
dependent. The most desirable solution would be to introduce the cutoff in some 
gauge independent way. These problems have been addressed by many authors 
and recently, significant progress has been made [46]. 

5.3 Alternative methods 

A variety of other techniques with which we can deal with the renormalization 
group have been developed by a host of other authors. Some of the correspond- 
ing numerical estimates are compared with those due to the Local Potential 
Approximation in section 2.11. Although a detailed analysis of these methods 
is beyond the scope of this thesis, for completeness we include some references 
to the e-expansion [1,16,42], lattice methods [47-49], perturbation theory [50-53] 
and large N expansions [1,42,43,54]. 
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